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ABSTRACT 


This thesis examines the quantitative difference between multi¬ 
ply and singly scattered infrared terrestrial radiation in an 
aerosol-laden planetary boundary layer (PBL). The total net 
(scattering plus emission) flux densities are calculated at 11 

levels i for each scattering process over a spectral range of 350 
-1 -1 -1 

cm to 2000 cm at 50 cm intervals .^The model PBL is a sum¬ 
mertime mid-latitude urban structure from LOWTRAN 7. The study 
prototype incorporates an inversion located at 2.00 km above a 
surface whose elevation is at sea level with a visual range of 
4.0 km. The visual range and extinction coefficients decrease 
with height to minimums located just below the inversion to 
simulate a differentially mixed aerosol-laden PBL. 

MUrfp * * S< « i ■ ’ ■ * •* »*•"'' (>»,*' • ■ y ' 1 - 

' Results show that a pronounced total downward directed net 

flux develops in the presence of aerosol particles in the water 
vapor absorption band between 1200 cm ^ to 2000 cm This 
negative total flux is 44 times larger then the singly scat¬ 
tered flux densities. It appears that a synergistic relation¬ 
ship between absorption/emission and multiple scattering among 


the radiatively active gases and the aerosols is established. 


In the water vapor window region of the infrared spectrum 
(830 cm * to 1250 cm *), thermal emission from the surface is 
able to penetrate to the top of the urban PBL. Here the differ¬ 
ence between multiple and single scattering is reduced but mul- 
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tiple scattering is still numerically superior ranging from 
0.02 to 1.02 times the single scattered flux densities. 
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This thesis examines the quantitative difference between multi¬ 
ply and singly scattered infrared terrestrial radiation in an 
aerosol-laden planetary boundary layer (PBL). The total net 
(scattering plus emission) flux densities are calculated at 11 
levels for each scattering process over a spectral range of 350 
cm ^ to 2000 cm ^ at 50 cm ^ intervals. The model PBL is a sum¬ 
mertime mid-latitude urban structure from LOWTRAN 7. The study 
prototype incorporates an inversion located at 2.00 km above a 
surface whose elevation is at sea level with a visual range of 
4.0 km. The visual range and extinction coefficients decrease 
with height to minimums located just below the inversion to 
simulate a differentially mixed aerosol-laden PBL. 

Results show that a pronounced total downward directed net 
flux develops in the presence of aerosol particles in the water 
vapor absorption band between 1200 cm 1 to 2000 cm *. This 
negative total flux is 44 times larger then the singly scat¬ 
tered flux densities. It appears that a synergistic relation¬ 
ship between absorption/emission and multiple scattering among 
the radiatively active gases and the aerosols is established. 

In the water vapor window region of the infrared spectrum 
(830 cm * to 1250 cm *), thermal emission from the surface is 
able to penetrate to the top of the urban PBL. Here the differ¬ 
ence between multiple and single scattering is reduced but mul- 
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DIGEST 


This thesis examines the quantitative difference between multi¬ 
ply and singly scattered infrared terrestrial radiation in an 
aerosol-laden planetary boundary layer (PBL). The total net 
(scattering plus emission) flux densities are calculated at 11 
levels for each scattering process over a spectral range of 350 
cm ^ to 2000 cm * at 50 cm * intervals. The model PBL is a sum¬ 
mertime mid-latitude urban structure from LOWTRAN 7. The study 
prototype incorporates an inversion located at 2.00 km above a 
surface whose elevation is at sea level with a visual range of 
4.0 km. The visual range and extinction coefficients decrease 
with height to minimum? located just below the inversion to 
simulate a differentially mixed aerosol-laden PBL. 

Results show that a pronounced total downward directed net 
flux develops in the presence of aerosol particles in the water 
vapor absorption band between 1200 cm * to 2000 cm This 
negative total flux is 44 times larger then the singly scat¬ 
tered flux densities. It appears that a synergistic relation¬ 
ship between absorption/eraission and multiple scattering among 
the radiatively active gases and the aerosols is established. 

In the water vapor window region of the infrared spectrum 
(830 cm * to 1250 cm ^), thermal emission from the surface is 
able to penetrate to the top of the urban PBL. Here the differ¬ 
ence between multiple and single scattering is reduced but raul- 
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tiple scattering is still numerically superior ranging from 
0.02 to 1.02 times the single scattered flux densities. 

It is concluded that for the total terrestrial infrared 
domain multiple scattering above and beyond single scattering 
must be considered. Furthermore, multiple scattering in that 
domain leads to a coupling between the redirection of photons 
with increased gaseous and particulate absorption and emission 
at the respective layer temperature. 
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1.0 


INTRODUCTION 


1.1 Overview 

Historically, the development of modeling radiative 
transfer in the atmosphere can be characterized as a gradual 
increase in complexity. Most early endeavors were concerned 
exclusively with short wave radiation and the phenomenon of 
scattering by molecules. The field began to expand in tandem 
with quantum mechanics and scattering theory. These branches of 
physics allowed investigators to include the effects of both 
scattering and absorption by a multitude of larger particles 
found in the atmosphere. But, even this advance kept the 
interaction of electromagnetic radiation with particles a 
purely "local" phenomenon concerned only with the individual 
particles reeraissions and single scattering interplay. Stymied 
by such constraints the true nature of aerosol scattering 
within the planetary atmosphere can not be fully appreciated. 
In order to more completely address the question of aerosol 
scattering and absorption, investigations must include the 
effects of multiple scattering and expand the electromagnetic 
spectrum to include the thermal planetary domain. 
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1.2 Definition of the Research Problem 

I propose to explore the question of how quantitatively 
significant is multiple scattering within the thermal long wave 
domain above and beyond single scattering in an aerosol laden 
boundary layer. Several cases will be investigated to determine 
the significance, if any, of the presence of particulate aero¬ 
sol matter for infrared multiple scattering. Physically, multi¬ 
ple scattering on aerosol particles increases the path length 
that a photon will travel, thus increasing the probability of 
the photon interacting with all the constituents in the layers 
under consideration. A comparative examination of single 
scattering versus multiple scattering in several scenarios 
should allow insight to be gained into the importance of each 
condition as it pertains to a single control case. In addition, 
the inclusion of the intrinsic infrared radiative transfer from 
the surface boundaries and atmospheric layers, previously not 
included in most studies, will be more realistic when surveying 
the upward and downward fluxes. 

1.3 Research Objectives 

My topical research objective is to compare, by juxtaposi¬ 
tion, the quantitative difference between multiply (higher 
order) and singly (first order) scattered infrared terrestrial 
flux densities in an aerosol-laden boundary layer. This will 
render an assessment about the significance, if any, of aerosol 
scattering in the longwave domain near the earth's surface. 
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One methodological research objective is related to my 
primary goal. This objective is the testing (verification and 
validation) and implementation of a computer code capable of 
simulating multiple scattering. To reach this objective, a com¬ 
plete restructuring of the delta-Eddington code, as presented 
by Wiscombe (1977) to solve for solar flux densities, to a code 
which is capable of producing multiply scattered infrared 


fluxes will be undertaken. 








2.0 REVIEW OF PERTINENT LITERATURE 


2.1 Infrared Radiative Effects of Aerosols 

There are important reasons for pursuing aerosol radiation 
studies. Some of these are: 1) reduced insolation to the tropo¬ 
sphere, 2) the addition of anthropogenic soot, with its high 
absorption capability, to the atmosphere on a truly global 
scale, 3) global climate models lack any aerosol influence in 
their radiation parameterizations, and 4) aerosol pollution of 
all sorts substantially alters the urban boundary layer through 
radiative-dynamic interactions. Excellent reviews by McCartney 
(1976) and Wiscombe (1983) chronicle the study of the radiative 
effects of aerosols. However, it is the last reason which is 
the impetus for investigations of the infrared effect of aero¬ 
sols within planetary boundary layers (PBL). Viskanta and Weir- 
ich (1979) drew the following conclusions about the effects of 
aerosols in the PBL. 1) Radiatively participating aerosols in 
the urban environment generally give rise to cooling of the 
ground and warming of the PBL. Whether there is a net cooling 
or warming of the earth-PBL system depends upon the surface 
albedo, the radiation absorption, scattering characteristics of 
the aerosol, aerosol loading, and the vertical distribution and 
zenith angle (latitude, season and time of day). 2) Radiatively 
active gaseous pollutants alter the local energy balance and 
through it modify the vertical temperature structure, 
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stability, and heat fluxes at the ground. 

In the early 1970's environmental concerns prompted an 
interest in the effects of gaseous and particulate pollutants 
on the radiation budget of the earth-atmosphere system on both 
global and regional scales. The regional studies have generally 
been concerned with modeling polluted urban boundary layers, 
and considering the effect of increased aerosol loading on the 
net radiative flux in the PBL. The scattering and absorption 
of visible (shortwave) radiant energy has been treated quite 
accurately by many investigators to include Welch and Zdun- 
kowski (1976), Viskanta et al. (1977), Coakley and Cess for 
climate models (1982), and Dave (1981). 

Treatment of infrared radiative transfer in the PBL, how¬ 
ever, has been considerably less comprehensive and complete. 
Wang and Domoto (1974) were one of the first investigators to 
employ a modified two-stream approximation to compute the 
transfer of infrared radiation in a finite inhomogeneous, tur¬ 
bid atmosphere. Ackerman, Liou, and Leovy (1975) used a four 
stream, multi-layered radiative transfer model to treat the 
transfer of infrared radiation in the atmosphere containing 
both scatterers and absorbers. They concluded that aerosols 
produce insignificant influences on temperature structures in 
bands which exhibit strong gaseous absorption, but their 
effects are of considerable importance in the window region. 
They also found for a typical case with a low lying inversion 


at about 1.0 km and a uniform aerosol concentration below that 
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cooling in this layer can be on the order of 5“K or 10 a K per 
day and that the aerosol effects in the infrared may influence 
the strength and positions of inversions. 

2.2 Infrared Multiple Scattering Methodologies 

Radiative transfer theory provides the mathematical 
description of the interaction between both external (incident 
solar) and internal (emitted thermal) sources of electromag¬ 
netic radiation and the radiatively active atmospheric consti¬ 
tuents of gases and particulates. The solution to the radia-r 
tive transfer equation for a scattering atmosphere (Chan¬ 
drasekhar, 1960; Goody, 1964; Liou, 1980) generally require a 
numerical' solution because of the integro-differential equa¬ 
tions involved and that the "exact solutions" are very limited 
in scope of applications. So much computer time would be 
required to solve these exact solutions that it would be 
impractical to use the methods routinely. However, the "exact 
solutions" they could be used as a check for approximate numer¬ 
ical methods. Survey papers by Lenoble (1977) and Hansen et al 
(1974) discuss the various multiple scattering routines both 
exact and approximate. However, explicit mention of a routine 
specifically addressing infrared radiative transfer is not 
present. The comment is only made that a thermal term may be 
added to the regular radiative transfer equation to account for 
the diffuse component of internally generated thermal radia¬ 


tion . 
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The most commonly used numerical solution tc the radiative 
transfer equation is the two-stream approximation. In general, 
the two-stream approach is a method of representing the diffuse 
radiation field in two streams in one-dimensional radiation 
transfer problems. The approximation proposed by Schuster 
(1905) and Schwarzschild (1906) is to assume that the intensity 
in the positive direction is isotropic and that in the negative 
direction the intensity has a different value but is also iso¬ 
tropic. Meador and Weavor (1979) provide a compendium of two- 
stream methods and observed that all existing two-stream 
approximations can be represented by identical forms of coupled 
differential equations if the radiation intensity is replaced 
by integrals of the intensity over hemispheres. The two-stream 
methods avoid the complex and lengthy computer procedures 
necessary for numerical solutions of the differential equations 
while yielding closed-form analytical results that are rela¬ 
tively easy to interpret. 

Shettle and Weinmann (1970) used an Eddington approxima¬ 
tion (a form of the two-stream method) to calculate the diffuse 
solar flux densities in models of homogeneous turbid layers. 
In their paper they suggested that the infrared domain may be 
breached by adding to the radiative transfer equation an 
appropriate term to represent the thermal emission from layers 
under scrutiny. The same information is provided by Liou 
(1980). However, neither paper derives explicitly the solution 
to such an addition to the radiative transfer equation. 
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Joseph and Wiscombe (1976) identified a problem with the 
Eddington approximation in the presence of anisotropic scatter¬ 
ing particles. Typically, the scattering of light in planetary 
atmospheres is in most cases characterized by a significantly 
anisotropic phase function (scattering diagram) since the pres¬ 
ence of even a small number of haze or cloud particles results 
in a forward elongation of the scattered radiation field inten¬ 
sity. They overcame this problem using a delta function to 
effectively truncate the forward diffraction peak produced by 
the asymmetric scattering pattern of aerosol particles. While 
this method expanded the Eddington technique to anisotropic 
scattering conditions, it was restrictive in the sense that it 
could only be applied to homogeneous atmospheres with the same 
optical properties. Wiscombe (1977) expanded the delta- 
Eddington method to inhomogeneous layers by adopting the 
parameterized values of Joseph and then solving the resulting 
differential equations by a matrix method solution. This was 
possible because he made the assumption that flux continuity is 
enforced at the layer interfaces. Wiscombe's methodology 
yielded very accurate diffuse solar fluxes for inhomogeneous 
turbid atmospheres for a broad range of optical depths. This 
condition would have to be met if the flux densities from an 
urban EBL were to be compared for single and multiple scatter¬ 
ing since the optical properties within the PBL change rapidly 
at various elevations. The simplicity of Wiscombe's approach 
and the accurate results of the method lead naturely to its 
choice as the model to convert to the thermal domain. 







3.0 BASIC THEORY AND DATA BASE 

Radiative transfer theory for aerosol scattering absorp¬ 
tion is based on the Mie theory for particles whose radius is 
greater than 0.03 times the wavelength of the irradiating 
electro-magnetic wave energy. Assuming time independence, 
elastic scattering (no conversion from one frequency to another 
in the range of observation), no internal sources (e.g., ther¬ 
mal emission from an embedded target, like a cloud), and suffi¬ 
ciently rare media that each particle is separated to the 
extent that no interparticle shadowing occurs, the radiative 
transfer equation for the diffuse radiance I(t,p,9), as lucidly 
presented by Shettle and Weinman (1970), is: 

U ^ (t,p,$) - -I(t,p,<j>) 

- 2ti *1 

+ I I P(p>$;m' )I(t,M 1 »' )dp' d<j>' (3.1) 

o -1 

_ z 

« _ , v P o 

+ r F o ) e 

4 O 0 0 

where the second term is the diffuse component i.e., multiply 
scattered radiance (self-illumination) and the third term 
represents the collimated direct solar beam duly attenuated to 
optical depth (t) by extinction. 

The following notations and Figure 3.1, Scattering 
Geometry, will expound upon the elements making up the terms of 
Eq.(3.1). 


9 
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Referring to Figure 3.1: 

|i - cos9 

9 - zenith angle 

$ - azimuth angle 

© = angle between incident and scattered radiances 

z - altitude 

(z) » Extinction coefficient implicitly present in the opti¬ 
cal depth (t). It is a bulk property of the medium derived from 
Mie Theory, and it is composed of four other descriptive quan¬ 
tities dealing with the attenuation process. 

1. ) molecular absorption (k ) 

di 

2. ) molecular scattering 0^ s ) 

3. ) particle absorption (o ) 

a 

4. ) particle scattering (o ) 

s 


z - optical depth given by integrating the relevant vertical 


extinction coefficient profiles: 


J (k +k +o +o ) dz 
a s a s 
z 

u « albedo for single scattering defined as the ratio of the 
o 

amount of flux scattered to that scattered and absorbed 

(k +o ) / (k +k +o +o ). 
s s a s a s 






Figure 3.1 Scattering geometry showing the relationship 
between the zenith, azimuthal, and phase function angles (after 
Liou, 1980). 
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( I 

P(p,<j>,p ,4> ) «* The phase function which defines the fractional 

I I 

amount of light incident at (p , $ ) which is scattered in the 
direction (p,4>). 

7iF q - Solar irradiance perpendicular to the direction of 
incidence. 

Mie Theory provides a "complete" solution to the problem 
of scatter and absorption by a particle represented as an iso¬ 
lated isotropic sphere. The scattering properties depend on the 
size of the sphere, the wavelength of the incident flux, and 
the real and imaginary parts of the complex index of refrac¬ 
tion. The size values are related in the Mie factor (a) defined 
by 

a = 2irr/X (3.2) 

where (X) is the wavelength and (r) is the radius of the 
sphere. It is the relative size of the particle in terms of 
wavelength which is most important for determining the scatter¬ 
ing process. When the particle is far smaller than the 
wavelength, the scattering is called Rayleigh scattering. 
Scattering of this type varies inversely as the fourth power of 
the wavelength. The principal Rayleigh scatterers in the atmo¬ 
sphere are the molecules of the constituent gases (a < 0.03). 

Cloud and aerosol particles are a different matter. Atmos¬ 
pheric particles fall broadly into four categories. (1) Solid 
aerosols of mean size between 0.1 and 1 pin; they are not neces¬ 
sarily spherical and their refractive index may be highly vari- 
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able. (2) Haze water drops of mean size between 0.1 and 1pm; 
they are spherical and have known composite refractive indices. 
(3) Cloud water drops of mean size between 1 and 100 pm; again 
spherical and of a known refractive index. (4) Cloud ice par¬ 
ticles of mean size between 1 and 100 pm; irregular shape but 
known refractive index. Therefore, the relative size factor 
(a) becomes even more important since the scattering pattern in 
three dimensions changes dramatically as the radius of the par¬ 
ticle increases along with changes in the absorption qualities 
of the aerosol particle in question. 

3.1 Single vs Multiple Scattering 

Single scattering assumes that the particle is exposed 
only to the light of an incident beam. No account is taken of 
the fact that each particle in a scattering volume is exposed 
to and also scatters a small amount of radiation already scat¬ 
tered by other particles. The LOWTRAN models prior to the 
release of LOWTRAN 7 only included the effects of single 
scattering and then only as a loss mechanism whereby scattering 
was not considered as a source mechanism for the path radiance. 
Therefore, single scattering only depleted the number of pho¬ 
tons from a scattering volume not allowing the input from more 
distant sources. 

Multiple scattering, on the the other hand, includes 
higher orders of scattering when each particle is exposed to, 
and also scatters light which has already been scattered by 
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other particles. Obviously, multiple scattering is an important 
process for the transfer of radiant energy in the atmosphere, 
particularly, when aerosols are involved. Multiple scattering 
allows the possibility of a global intervention from particles 
in the far field. This multiplicity enhances the probability 
that a given photon will be absorbed and/or scattered. This is 
so because the effective path length is increased in the pres¬ 
ence of radiatively active gases and aerosols. 

3.1.1 Absorption-Emission 

Scattering is often accompanied by absorption. Both 
processes remove flux from a given beam of light. Scattering is 
explained'in terms of the wave theory of light and it produces 
no net change in the internal energy states of the molecules. 
Absorption, however, requires quantum theory for its explana¬ 
tion and does produce changes in energy states. Three forms of 
internal energy exist: (a) rotational, (b) vibrational, and (c) 
electronic states. The incident radiation is regarded as being 
quantized and as such only whole quanta can be accepted by the 
absorbing molecule. In absorbing a quantum of energy the 
molecule undergoes a transition from a lower to a higher state 
of one of the three internal energy forms. When absorption is 
significant for a particular aerosol, this fact is expressed in 
the imaginary portion of the complex index of refraction. The 
larger this number is the greater the amount of energy which 
can be taken out by absorption. 
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Absorption is only the first part of a cycle which is com¬ 
pleted by emission. As a consequence of molecular motions and 
collisions, molecules endlessly exchange internal for transla¬ 
tional energy, and vice versa. Molecules already excited to 
upper levels are "relaxed", while molecules at lower levels are 
"excited" to upper levels by collisions with both events 
operating on a statistical basis. Most of the upper levels are 
inherently unstable, however, and molecules occupying these 
levels undergo a transition to lower levels by emitting quanta 
of radiant energy. The above physical activities form the basis 
behind blackbody radiation theory. Emission is as spectrally 
selective as absorption and is a discontinuous process. The 
difference between the flux removed from an incident beam and 
the totally scattered flux must be attributed to absorption by 
the particle. The absorbed energy goes into heating the parti¬ 
cles which as described above leads to further emission. A 
cycle in which the medium surrounding the particles is heated 
leading to changes in the ambient environment due to the pres¬ 
ence of aerosols. 

3.2 Bulk Properties From Mie Theory 

In the field of light scattering and radiative transfer, 

it is customary to use a terra called cross-section, which is 

analogous to geometrical area, to denote the amount of energy 

removed from the beam by particles. In the case when the 

2 

cross-section refers to a particle, its units are in area (cm ) 
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as given in cgs units. Thus, the extinction cross-section, in 

units of area, is the sum of the scattering and absorption 

cross-sections. However, when the cross-section is in reference 

2 -1 

to a unit mass, its units are in area per mass (cm g ). The 

mass extinction cross section is therefore the sum of the mass 

absorption and mass scattering cross-sections. Furthermore, 

when the extinction cross-section is multiplied by the particle 

_3 

number density (cm ) or when the mass extinction cross-section 

-3 

is multiplied by the density (gem ) the quantity is referred 
to as extinction coefficient, which has units of length either 
(cm ) or as with LOWTRAN (km ^). In the field of infrared 
radiative transfer the mass absorption cross-section is simply 
referred to the absorption coefficient as cited by Liou (1980). 


Since the particles occupying a particular volume normally 
have a size distribution (such that the total number of parti- 

a> 

cles n » / n(r) dr where n(r) is the number of particles with 
0 

radii from (r) to (r+1), the mass or (volume) extinction coef¬ 
ficients and the phase function of a unit volume of the medium 
must be calculated from appropriate integrals over all the par¬ 
ticles. That is: 


B - I J(r) n(r) dr 
sc 0 
oo 

B - f K(r) n(r) dr 
6 0 
09 

P(9) - 1/S / n(r) J(r) P (r,9) dr 
0 


(3.3) 

(3.4) 

(3.5) 


where J(r), K(r) and P (r,8) are the scattering and extinction 








coefficients and the phase functions appropriate to the indivi¬ 
dual particles of radius (r) (Paltridge and Platt, 1976). 

3.2.1 Efficiency Factors and Refractive Index 

The single particle scattering and extinction coefficients 
J(r) and K(r) are normally given in terms of an expression for 
the corresponding efficiency factors Q and which are sim¬ 
ply the ratios of J(r) and K(r) to the geometric cross section 
2 

nr . Thus according to Paltridge and Platt (1976): 

a> 

Q - J(r)/nr 2 = 2/a 2 I (2n+l) (a a*+b b*) (3.6) 

sc , n n 

n» 1 

and 

oo 

Q = K(r)/rrr 2 =* 2/a 2 E (2n+l) Real (a +b ) (3.7) 

e , n n 

n=l 

where the a and b are complex quantities, taken from Mie 
n n 

theory, which depend on the size parameter (a) and the complex 
index of refraction (m), 

When absorption occurs the index of refraction becomes 

ra « m - ira. (3.8) 

r l 

with and m^ representing the real and imaginary parts of the 
complex index of refraction, and (i) the imaginary unit. The 
imaginary portion carries the value of the amount of absorption 
taking place when electro-magnetic radiation interacts with a 
material substance. These are usually laboratory measurements 
under controlled conditions. Figure 3.2 shows the scattering 
efficiency factor Q sc as a function of the size parameter (a) 
for m^ - 1.33 with several differing values of the imaginary 
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Figure 3.2 Efficiency factor for scattering, Q , as a function 
of the size parameter a - 2nr/X.. The refractive index in 
ra -1.33, with results shown for four values of m. (after Hansen 
and Travis, 1974). 1 
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portion. For nu - 0 ,and the condition of a perfect reflector, 

there is no absorption, so that Q - Q . Q in this case 

sc e sc 

shows a series of major maxima and minima with ripples. The 

major maxima and minima are due to the interference of light 

refracted and transmitted by the sphere, whereas the ripples 

arise from edge rays that are grazing and traveling along the 

sphere, spewing off energy in all directions. Q or Q 

sc e 

increases rapidly when the size parameter reaches about five 
and approaches an asymptotic value of two there after. This 
implies that a large particle removes from the incident beam 
exactly twice the amount of light that it can intercept 
geometrically. Physically, the removal of the incident light 
beam includes the diffracted component, which passes by the 
particle, plus the light scattered by reflection and refraction 
within the particle. Both the ripples and the major maxima and 
minima damp out as absorption within the particle increases. 

3.2.2 Phase Function and Single Scattering Albedo 

The phase function for scattering expresses in a formal 
manner the angular dependence of scattering on the "viewing" 
angle. It is defined by van de Hulst (1957) "as the ratio of 
the energy scattered per unit solid angle in this (a given) 
direction to the average energy scattered per unit solid angle 
in all directions." This is symbolically represented either by 

P(9) - P(p,<t>;p',$') (3 * 9) 

where as stated earlier the 8 is the angle between the incident 
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and scattered fluxes (scattering angle). 

Related to the phase function is the asymmetry factor (g 

or <cos9>). This quantity physically represents the difference 

between the flux densities in the forward direction and the 

backward direction arising from the scattering by a volume (or 

particle) when the incident flux is normalized to unity. This 

forward minus backward over the total scatter ratio is defined 

as the integral over all solid angles of the phase function. 

That is the first moment expression: 

♦1 

<cos9> »■ / P(0)cos(0)d(cos9). (3.10) 

-1 

It is the minimal information required about the phase function 
which allows some accuracy in the approximate solutions of a 
transfer problem. It is equal to zero for isotropic scattering 
(i.e., the same intensity over both hemispheres), and to +1 and 
-1 for complete scattering into the forward and backward direc¬ 
tions, respectively. 

The asymmetry parameter is particularly useful when multi¬ 
ple scattering effects are significant because the details of 
the scattering phase function tend to be smoothed out. This is 
when the asymmetry parameter is used to characterize the angu¬ 
lar distribution of the radiation field since approximate 
methods like the Eddington and two-stream assume isotropy of 
the radiance in the forward and backward directions to model 
the azimuthally averaged phase functions pertinent to these 


methods. 
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3.3 Two-stream and Eddington Multiple Scattering Approxima¬ 
tions 


The Eddington and its close relative the delta-Eddington 
approximations are both members of the two-stream family of 
approximations for multiple scattering. My reasons for choosing 
the Eddington method was based on the following favorable 
attributes: (1) It provides some qualitative insight into the 
mechanisms whereby radiation is transferred within turbid 
atmospheres of moderate or large optical depths. (2) It allows 
rapid calculation of flux densities Ft,?.* within realistic 
inhomogeneous atmospheres with an accuracy of several percent 
compared to more exact methods like the discrete-ordinate 
method of Liou (1973). (3) The method's simple approximation to 
a complex physical process allows numerification by standard 
matrix inversion techniques as suggested by Wiscombe (1977). 

For all the two-stream approximations, three assumptions 
are made which help simplify the formulation of the governing 
equations. First, the transfer of radiation will occur in a 
horizontally plane-parallel atmosphere. Second, the radiation 
emission will be isotropic. Third, the upward and downward 
radiances, I + and I , at each level are replaced by the herais- 
pherically integrated average radiance given by the zeroth 
order expression 

I* - I I p dp. (3.11) 

0 ° 

This allows a complex scattering pattern to be represented by 
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an average intensity and when integrated over azimuth an iso¬ 
tropic quality. Taken together they allow a method whereby the 
radiation field can be specified in terms of only the upward 
and downward components of the flux density, Ft and F+, respec¬ 
tively. The irradiances are then directly related to the radi¬ 
ance by: 

Ft - rr I + 

(3.12) 

Fi - it i”. 

The Eddington approach assumes that the radiances can be 

given by a Taylor's series expansion, vis. 

I(t,p) - I(z) + I.(t)u (3 - 13) 

o 1 

where all higher order terms are neglected. Figure 3.3 shows 
visually ■ a comparison between the two methods together with 
their respective expressions for the flux densities. Each 
method contains its own phase function in terms of the asym¬ 
metry factor <cos9>. Thus for the two-stream case (Fig 3.3b) 

P(0) - l±<cos9>, (3.14) 

and for the Eddington approximation (Fig 3.3c) 

P(9) - !+3<cos9>cos(9). (3.15) 

As stated by Paltridge and Platt (1976), "Despite the fact 
that the two-stream and Eddington approximations are limited in 
applications to cases where the scatter is not highly anisotro¬ 
pic, they are highly valuable because of the simplicity of 
their analytical solution. Further, they retain a measure of 
directness in physical conception which allows easy and unambi¬ 
guous formulation of the boundary conditions in specific 









Two stream 



Ft = TT It Ft=TM+ I = I o + l,cos0 

F j, = TT I F TT I _ F*7T(I 0 ± §• I ,) 



Edd ington 


(a) (b) (c) 

Figure 3.3 Two-stream and Eddington approximations. (a) 
equivalent, (b) and (c) are conceptually equivalent (after Pal- 
tridge and Platt, 1976). 
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cases." 

3.4 Pertinent Data Base 

The primary optical values needed to develop a multiple 
scattering two-stream model with a minimum of inaccuracy are 
the optical depths (t), single scattering albedos (w), and the 
asymmetry factors (g) for each relevant layer. As defined above 
these values explicitly depend upon the efficiency factors out¬ 
lined previously. Few aspects of scattering have greater impor¬ 
tance or have attracted greater attention than atmospheric 
scattering and extinction and the associated matter of visibil¬ 
ity. Table 3.1 gives a listing of some of the published effi¬ 
ciency factors applicable to atmospheric particles. Similar 
sources exist for molecular scattering. Specifically, Elterman 
(1968,1970) tabulated values of the volume coefficients for 
many altitudes from sea-level to 50 km. My primary source for 
the optical properties came from a compilation by Shettle and 
Fenn (1979) which centered on lower atmospheric aerosol models. 
The aerosol properties were based upon experimental measure¬ 
ments that were made during and prior to the mid-1960's. The 
optical properties of the models are given for a number of 
wavelengths between 0.2 and 40.0 pra and for different relative 
humidities ranging from o to 99%. Specific to my topical objec¬ 
tive, the investigation of multiple versus single aerosol 
scattering in the terrestrial infrared domain, I chose a sum¬ 
mertime mid-latitude urban model for investigation. 









Reference 

Refractive 

index* 

Size 

parameter 

Hoi! (1948) 

1.33 

4.8. 5.4. 6.0. 

6.6 (24 values) 8.0 

Lowan (1949) 

1.33, 1.44. 1.55, 

2.00. 1.50 

0.5 (15 values) 6.0 

0.5 (24 values) 12.0 

Houghton and 

Chalker (1949) 

1.33 

7 (33 values) 24 

Gumprecht and 

1.2, 1.33, 1.4, 

1 (many values) 8. 

Sliepcevich (1951b) 

1.44, 1.5, 1.6 

10(5)100(10)200 

(50)400 

Johnson et al. (1954) 

Complex values, 

3.6 < A < 13.5 fjLtn 

1 (many values) 25 

Pangonis et al. (1957) 

1.05(0.05)1.30 

0.2(0.2)7.0(1)15 

Havard (1960) 

Complex values, 

3.6 < A < 13.5 n .m 

For particle r = 

1, 2. 4, 9, 12. 15 /zm 

Deirmendjian (1963b) 

Complex values. 

0.5(0.517.0, 10.0, 

8.15 < A < 16.6 nm 

15.0 

Penndorf (1956, 1957c) 

1.33, 1.40. 1.44, 

1.486, 1.50 

0.1(0.1)30.0 

Giese et al. (1961) 

1.5 

0.2(0.2)159 

McCormick (1967) 

1.5 

0.01(0.01)2.0(0.1) 

181.9 

Irvine and 

Complex values, 

For particle r = 

Pollack (1968) 

0.7 < A < 200 fim 

0.3, 1,3. 10 fim 

Zelmanovich and 

Shifrin (1968) 

Various complex 
values 

0.5 (many values) 100 


Table 3.1 Published tabulations of efficiency factors applica¬ 
ble to atmospheric particles (after McCartney,1976). 
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The attenuation coefficients for the different aerosol models 
(urban, rural, maritime, and tropospheric) are given in a 
tabular form and allow extrapolation for the desired wavenumber 
domain. An example of such a table is presented below in Table 


3-2. 
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4.0 


PROCEDURAL RESEARCH METHODOLOGIES 


4.1 Overview 

The topical objective of comparing aerosol single scatter¬ 
ing versus multiple scattering in the terrestrial long 
wavelength domain requires the development of a suitable algo¬ 
rithmic procedure to calculate the desired fluxes. Shettle and 
Weinman (1970) developed an Eddington formulation which dealt 
exclusively with solar fluxes. Their method was extended to the 
infrared regime by manipulation of the basic radiative transfer 
equation. A computer code analogous manner to Wiscombe's (1977) 
was constructed based on the delta- Eddington approximation of 
Joseph et al. ( 1976). 

This numerical model produces flux densities. However, 
validation procedures and comparative verification indicate 
that my development is not suitable. In search for a modified 
approximation, I turned to the recently released LOWTRAN 7 
which does contain a multiple scattering routine for the 
longwave domain. The insight gained from the first endevour, 
allows me to see the superiority of the LOWTRAN 7 two-stream 
and adding method. It is to be my primary tool into the inves¬ 
tigation of the scientific problem. 
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4.2 Derivation of Eddington Modification 

The development of a computer code suitable for determin¬ 
ing multiply scattered thermal radiation begins with a sugges¬ 
tion by Liou (1980). He indicates that for the transfer of 
terrestrial infrared radiation in scattering atmospheres, which 
are in local thermodynamic equilibrium, it is possible to 
replace the last term in Eq. (3.1) by a term representing a 
thermal emission contribution, vis., 

(1-5 v )[B v (T(t))] O-D 

where w is the single scattering albedo and B is the Planck 

V V 

function of temperature (T). The subscript v, dropped later, 
signifies that a spectral density function is explicitly under¬ 
stood. The temperature (T) is also dependent upon the optical 
depth (t) which in this development is taken to be the vertical 
coordinate for the used reference system. Using Eq. (4.1) in 
Eq. (3.1) yields: 

p ~ (t,p,$) = -I(z,v,4>) 

- 2n ♦ 1 

+ I I P (u, 4>; M ' ,<j>' )I(t,p' )dp'd<f>' (4.2) 

o -1 

+ (1-5 v )[B v (T(t))] 

where the wavenumber dependence is understood. Term one on the 
right, I(r,p,<j)), is the spectral radiance of the diffuse radia¬ 
tion field. Term two represents the multiply scattered diffuse 
component, and terra three is the thermal emission contribution 
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from optical depth t. Eq. (4.2) becomes the starting point for 
a procedure outlined by Shettle and Weinman (1970) to construct 
an Eddington approximation leading to two ordinary differential 
equations which when solved analytically produce simple 
exponential solutions for the flux densities. 

Eddington's approximation assumes the radiance can be 
given by 

I(r,p) = I (t) + pi . (4.3) 

o 1 

If the phase function is expanded as a series of associated 
Legendre functions, .in the manner of Chandrasekhar (1960), and 
that Eq. (4.3) represents the diffuse radiance, all terms of 
orders greater than one will vanish when Eq. (4.1) is 
integrated over (p) and (<j>). The phase function then may be 
approximated by 

P(cos(9)) = l+w^t) cos(0) (4.4) 

where 9 is the angle between the incident and scattered radi¬ 
ances (scattering angle). The integral in Eq. (4.1) then 
becomes 

2tt 1 

I i P(p,4>;p , ,<|>')I(p , ,<t> , )dp'd<t> < = 4ir(I +pgl ) ( 4 . 5 ) 

0-1 u x 

where 

w. *1 

g - T- - 1/2 / P(0) cose dcos(9). (4.6) 

-1 

is the asymmetry factor. Using Eq. (4.6) in Eq. (4.1) after 

integrating over and dividing by 2n produces 
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The particular solution will include a Planck function which in 
its unabridged wavenumber form is given by Walker (i960) as 

B v (T)-c 1 v 3 (exp(-c 2 v/T)-l). (4.13) 

2 

Here c^ is (2hc ) where (h) is Planck's constant and (c) is the 
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speed of light in the vacuum. This term in cgs units equals 
-12 2 -1 

1.1909x10 watts cm sr . is made up of Planck's constant, 
the speed of light and Boltzmann's constant combined algebrai¬ 
cally into (hc/k) to yield in cgs units 1.438 cm K®. However, 
numerically, this form proves to be cumbersome in the deriva¬ 
tion of the differential equations. I decided to follow Wallace 
and Hobbs (1976) and use a truncated Planck expression given 
by: 

B v (T(t>) -c 1 v 3 exp(-c 2 v/T). (4.14) 
This is done because throughout most of the wavenumber domain 
the exponential term in the untruncated version is normally 
much larger than unity. Only at the smaller wavenumbers (> 850 
cm *■) does the full expression become necessary. Testing indi¬ 
cates that a 4.0% error is incurred over the temperature range 
used in the models. Therefore, Eq. (4.11) becomes 


d 2 I /dt 2 -K I * -K c.v 3 (exp(-c,v/T(t)). (4.15) 

ooooi i 

cx 

The RHS of Eq. (4.15) is of the form (K e ) where (K) and (c) 
are constants and (t) is the independent variable. This implies 

CZ 

that the particular solution is simply some multiple of e 

Assuming that the particular solution I is equal to 

—c v/T(t) 

F e -I (4.16) 

and 


K. - -K c.v 
1 o 1 

with substitution into Eq. (4.15) yields 


(4.17) 


-c v/T(t) -c v/T(t) „ -c v/T(r) 

K.e - (F e ) - K F e 

1 o 


(4.18) 
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Implementing the 2nd derivative and solving for F along with 
substitution for produces 


3 3 4 

K c.v T K c.vT 

o 1 _ o 1 

2c 2 (4T/At) 2 c^AT/Ac) 2 


(4.19) 


which has units of radiance, watts cm , when multiplied by 

“ c 2 v/T 

e 


The AT/Ac elements in Eq. (4.19) represent the differential 
dT/dc which is implicit in the derivative of Eq. (4.18). No 
explicit relationship c^n be easily obtained for dT/dt because 
of the non-linear dependency of the extinction coefficient on 
the number density, gas pressure, and temperature. Therefore, 
a linear approximation is made with the input parameters AT and 
Ac determined from the LOWTRAN output. The change in tempera¬ 
ture (AT will come from arbitrary established levels and the 
changes in optical depth Ac will come from the wavenumber 
dependent mean transmittances (Tr) given by: 

Ac - -In (Tr). (4.20) 


The general solution for Eq. (4.15) is 


I 

o 



(4.21) 


where and 


are unknown arbitrary constants. 
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Similarly, a solution for 1^ yields 


A e 


A t 

o 


+ B e 


-A z 
o 


+ D e 


-c 2 v/T 


(4.22) 


where 

3(l-u))(l-ug)c 1 v 3 T 
° " 2c 2 (AT/At) 

3(l-u)(l-ug)c 1 v 2 T 2 
+ c 2 (AT/Ar) 

4-2 

+ e lV T 

(1-ug)(AT/At) 

which also has units of radiance when multiplied by e 


(4.23) 


c 2 v/T 


4.2.1 Computer Code Adaptation of Eddington Approximation 

Joseph, Weinman, and Wiscombe (1976) discovered the excel¬ 
lent accuracy of the delta-Eddington approximation for solar 
radiative flux (which calculates flux directly) for all phase 
functions, no matter how asymmetric, and for all optical depths 
and single scattering albedos. However, that work only inves¬ 
tigates homogeneous layers. Wiscombe (1977) developed an inho¬ 
mogeneous multi-layer version. I followed Wiscombe's procedure 
to develop a code which could generate fluxes for the thermal 
domain based upon the Eddington approximation of Shettle and 
Weinman (1970). The solution for solar flux in the delta- 
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Eddington method is equivalent to the Eddington approximation 
without the transformed parameters g',u', and t' of Joesph et 
al (1976). Thus, I feel justified in attempting to expand 
Wiscombe's numerical procedure to longer wavelengths. 


Accordingly, the diffuse irradiances may be computed from 

I and I. from 
o 1 

F* - n [I (t) ±2/3 I,(t)], k 1 * ,2U) 

o 1 

or in Wiscombe's notation 

G - n I , H - 2/3 tt I. (4.25) 

o 1 

from which 


Ft - tt (I -2/3 I,) - G - H 
o l 

F± » ir(I +2/3 I, ) =* G + H 
o 1 


(4.26) 


If the atmosphere is subdivided in such a way that the optical 
properties: optical, depth, single scattering albedo, and asym¬ 

metry factor all are constant within any layer then the com¬ 
plete layer structure can be modeled as a concatenation of N-l 
homogeneous sublayers as depicted in Figure 4.1. Here level 
(i+1) is topped by layer (i) with the vertical coordinate being 

"k 

the optical depth whose value is 0 at t(0) and t (N) at level 
N. By imposing top and bottom boundary conditions and requir¬ 


ing flux continuity across interior levels (interfaces) a sys¬ 


tem of simultaneous differential equations can be developed 
represented at layer (i) as 





s 


level 1 


level 2 


level 3 



w 1 »g 1 » (At) x 


> §2 * (At) £ 



= 0 


T 


2 



level N-l 

level N 


Lu N-l ,g N-l ,(AT ' ) N-l 

7T7TTTTT7 


T 


N-l 



Figure 4.1 Schematic presentation of Eddington modification 
atmosphere which is a concatenation of N homogeneous layers 
with layer values of single scattering albedo (w.), asymmetry 
factor (g.) and layer optical depth Lz. z^ is the total 
optical depth at level N with S being the incident solar flux 
at N-l. 
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G i - « [x 2 i _ 1 exp(X. i t i ) 
+ x 2i exp(-k.t.) 
+ a exp(-c 2 v/T)] 


(4.27) 


and 


H. - — 2/3tt [x,,. . exp(\.t.) 
1 2i-l ii 

+ x 2 . exp(-X i r i ) 

+P i exp(-c 2 v/T)] 


(4.28) 


where a. and fi. are from Eqs. (4.19 and 4.23) and X.. is >/k 
11 M 1 o 

from Eq.(4.10). The unknown constants for layer i are written 
in the notation 

*21-1 ' A i' *21 ' B i 

in order that the usual matrix/vector notation (C x « D) for a 
linear system of equations might subsequently be applied. The 
construction of the matrix system of equations rests upon the 
assumption of flux continuity at the interface between layers 
and solving for the unknown constants. Appropriately, the inte¬ 
rior flux continuity conditions are represented by 


V2 ( vi> - Vi'vi' ■ 0 

H n-2<Vi> - ■ 0 


(4.29) 


which form a couple for each layer, where t„....-c , are the 

L n-1 

level values of optical depth from the top to level N. Each set 
must be solved to find the unknown constants x 2^_2 and x 2i' T ^ e 
flux continuity conditions Eqs. (4.29) are expressed in terms 
of G and H rather than the flux density terms (Ft- G-H) and 
(Fi-G+H). This is because, by doing so, the coefficients have a 
simpler form, and numerically ill-conditioning is avoided. 
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Attempts to solve Eqs. (4.29) directly would lead to 

insurmountable ill-conditioning for the layers, no matter how 

they were subdivided. This is directly traceable to the 

(-\.t.. 

exponentials e and e in Eqs. (4.27 and 4.28) 

which become very large and very small for increasing optical 
depths (t^), respectively. Wiscombe observed that the flux den¬ 
sity formulas at each level required, not the bare quantities 

x. but rather 
i 


21-1 


X.. t. 
l I" 


i i*l 


2i-l 


(4.30) 


which come about because of the requirement for flux continuity 
across the interfaces. Substitution of Eqs. (4.30) into (4.27) 
and (4.28) produce the following interior flux continuity equa¬ 
tions : 


X.._, + x„- 


X„ . , 

2i* 1 


‘2i-l ' "2i exp(\. .At. ,) X 2i»2 eXp ^i-l At i*l) 

i*l i+l 


- a. +1 exp(-c 2 v/T. +1 )- a. exp(-c 2 v/T.) 


(4.31) 


and 


* « 2 i ♦ 1 a 

x„ . . + x„ . - -r--7-r - x„ . „exp(\. .At. .. 

2i-l 2i exp(X. .At. ,) 2i*2 r i*l i*l) 

i*l i*l 

- 3 i+1 exp(-c 2 v/T. +i )- 0. exp(-c 2 v/T i ) > 


(4.32) 


However, this substitution will prove to be fatally flawed 
since then the optical properties to be considered would only 
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come from the layer immediately adjacent to the level under 
consideration. Thus not allowing the other layers to leave 
their "imprint" on the layers under scrutiny and not attenuat¬ 
ing the transmitted emission. Eqs. (4.31) and (4.32), and as 
yet unspecified boundary conditions form the coefficient matrix 
C to solve for the unknown constants x 2 ^_^ and in C x * D 
where D represent the constant vectors of the RHS's of Eqs. 
(4.31) and (4.32). The requisite interior flux formulaes at 
each level (i+1) would be: 


F Vl “ 5/37T t x 2i-l. 


k 2i 


+ O i + cr) exp(-c 2 v/T i ) ] 


(4.33) 


and 


Fi - , 

l+l 




+ (a i -0 i ) exp(-c 2 v/T.)] 


(4.34) 


4.2.2 The Boundary Conditions 

The top boundary condition poised a unique problem since 
the assumption of local thermodynamic equilibrium (LIE) was 
made. LTE does not occur within the atmosphere till approxi¬ 
mately 45 kilometers based on considerations presented by Goody 
(1964). Therefore, the emissions from the layers above have to 
be adjusted so as not to project the full Planckian emittance. 
Choosing as does LOWTRAN the top of the atmosphere at 100 
kilometers the Planckian function had to be scaled to perform a 


smooth transition to levels where LTE sets in. A 









40 


proportionality constant will be multiplied by the Planckian 
emission from that level reducing the output. For example, at 
100 km, the constant would be As a boundary condition input, 
the summation of the levels would constitute a "direct" ther¬ 


mal (F ) input from layers which although not at LTE, still 
o 

produce an emission due to the presence of a tenuous atmo¬ 


sphere. The top boundary condition becomes 


5/3tt - 77 ^-—. ce*x„ . exp(k.c 0 ) 

expCX^t^ 2 i 1 2 

- F q - exp(-c 2 v/T )(a +3 l ). 


(4.35) 


The bottom boundary condition uses the common assumption 
of a Lambertian reflectance, i.e., the upward isotropic flux 
given by a constant surface albedo (A) multiplied by the down¬ 
ward flux, plus the emission from the underlying surface 
represented by: 

3 

c.v exp(-c„v/T ) 
l t s 

produces the upward flux density from the surface. But, for a 
complete description, the layers approximately 1.0 km above the 
surface boundary are also included as suggested by Pallmann 
(1988) based on tables found in Goody (1964). This is accom¬ 
plished by multipling the emittance from each layer within this 
range and duly attenuating the cumulative radiance by the opti¬ 
cal depth of each layer. This action then represents diffuse 
counter-radiation which can be significant in aerosol-laden 
environments. The formulation then leads to an analogous 
expression for the bottom boundary condition used by Wiscombe 
for the solar regime to include the reflected radiation 
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(reflux), vis., 

x 2n _ 3 t (l-A)ti-(l +A )2/3n) ] 

+ x 2n _ 2 [ (1 ~ a )tt-( 1 * a )2/3t t ] 

3 (4.36) 

- e c.v exp(-c.v/T ) 

1 is 

+ reflux ( 1 -e). 


4.2.3 Validation and Verification of the Eddington 
Thermal Modification 

Upon activating the algorithm and programming actual runs 
of Eqs. (4.31) and (4.32) as submitted to the original penta- 
diagonal L-U decomposition matrix inversion method prove to be 
ill-conditioned for the generation of a non-singular result. To 
solve this problem I adopted a standard singular value decompo¬ 
sition method from-Press et al. (1986), specifically designed 
for systems which approach the condition of singularity. Ill- 
conditioning of the system is avoided but the values for the 
fluxes are unacceptably too large when compared to the LOWTRAN 
6 single scattering routine and the direct thermal emission 
based on the full and abridged Planck funtion. Both methods 
indicate that ray results were anywhere from four to five orders 
of magnitude too large. 

I believe the reason for this is that my development of 
the Eddington differential equations failed to produce an ade¬ 
quate physical model of the real attenuation process taking 
place in optically thick environments. The adjacent layer opti¬ 
cal thickness, transmittance, and reflectance were not present 
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implicitly or explicitly in my formulation. The only input to 
the individual level flux values came from the layer immedi¬ 
ately below. Thus, Wiscorabe's original substitution for the x 
values to avoid ill-conditioning in the matrix proved not to 

A 

take into account the total optical depth (t ) at the level 
under scrutiny i.e., the fluxes became isolated values. The 
upward and downward flux values, not duly attenuated, become 
inflated and do not represent the mulitple scattering taking 
place. It is obvious that this method is not appropriate and a 
replacement must found. Due to time constraints, LOWTRAN 7, the 
latest version of the LOWTRAN family and-newly released, is 
chosen. 

4.3 LOWTRAN 7 

4.3.1 Overview 

Although the Eddington modification did not prove fruitful 
it allowed insights into the complex nature of multiply- 
scattered radiation fields. The approximation, suggested by 
Shettle and Weinmann proved insufficient to account for all the 
physical processes which occur in the atmosphere in the pres¬ 
ence of a heavily aerosol laden planetary boundary layer. 

LOWTRAN 7 is capable of handling the complex multiple 
scattering problem by using an adding method and a k- 
distribution function. The k-distribution function adequately 
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decouples the scattering of aerosols from the absorption by 
gases allowing a separation of the two phenomena. This step is 
necessary since LOWTRAN 7 is a low spectral resolution model 
which does not allow line by line integration of the molecular 
absorption bands. The adding method algebraically combines 
layers to include the processes of transmission and reflection 
to account for the flux densities emanating upward and downward 
from the layer interfaces. 

The following outline of the mechanics of LOWTRAN 7 ori¬ 
ginated from a technical report authored by Issacs et al. 
(1986). Although, this report was written for adaptation of 
LOWTRAN 6, the same formulations were used for both models. The 
other papers referred to in this section are used to clarify 
the applicable procedures and are duly cited. 

4.3.2 Stream Approxiraation/Multiple Scatter 
Source Function 

In the stream approximation, the multiple scattering con¬ 
tribution to the source function J 

o 

J (t.p,*) - u> o /4IT / P(0;Q')I(t,n-)dO' (4.37) 

Q' 

is approximated by assuming isotropic scattering radiances I + 
and I over upward (Q + ) and downward (Q ) hemispheres, respec¬ 


tively. 
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Substitution of the azimuthally averaged radiances 
(4.37) yields: 

J M _(r,ii,*) - co /4ir [I + (c) / P(fi,fi + )dO + 

Wo O 

♦ 

a 


into Eq. 


(4.38) 


+ l“(t) I P(Q,Q - )dn ]. 
Q~ 


This approximation is the same as that used for the Eddington 
method based on isotropy over the appropriate hemisphere. 
Integrating over the angular scattering functions for the 
resulting azimuthally averaged backscatter fractions, f$(p,g) 
(based on a function of zenith angle cosine and the asymmetry 
factor), and substituting the corresponding fluxes I±(t) = 
F±(t)/it results in the multiple scattering source funtion to 
be : 


J M -(c,iy,) - co /it [F±(t)(i-P) + F±(t) $]. -39) 

The quantity (5 is the backscatter fraction for isotropically 
incident radiation and is based on elliptical integrals involv¬ 
ing the zenith angle and asymmetry factor as presented graphi¬ 
cally in Fig (4.2). This method uses a Henyey-Greenstein (H-G) 
phase function and in doing so introduces a small amount of 
error (< 32). Howe er, this small error is quite acceptable 
compared to the advantage Hansen (1969) demonstrated when he 
showed that H-G function can be used to replace the more real¬ 
istic Mie phase function in multiple scattering calculations 


with no more than a few percent error. 






BACKSCATTERED FRACTION 



ASYMMETRY FACTOR, q 


Figure 4.2 Backscattered fraction B and B(p) for the Henyey- 
Greenstein phase function versus the asymmetry factor (g) for a 
range of values of p (Wiscombe and Grams, 1976). 
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The source function is integrated along the desired path to 
obtain the desired total radiance including, now, the approxi¬ 
mated multiple scattering contribution. Here the fluxes are 
taken as the layer mean quantities evaluated at a level halfway 
through the layer. The asymmetry factor, (g), is a measure of 
directional scattering and is available from tabulations as in 


Shettle and Fenn (1970). 
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4.3.3 Layer Fluxes 

Fluxes approximate the required radiances for evaluation 
of the multiply scattering source function. There are three 
types of fluxes. The first are the "local" fluxes (F ± )which are 
produced by each individual local layer. Next, come the "compo¬ 
site" fluxesCF^*) which are produced by the adding method when 

individual intrinsic layers are combined. Finally, the "total" 
2 ± 

fluxes(F ) are generated when the adding algorithm makes its 
upward and downward passes through the model atmosphere layer 
by layer. 

The local fluxes come from the flux parameterizations for 
each ' individual intrinsic layer. Included in the parameteriza¬ 
tion are the transmission (T) and reflection (R) functions 
based upon the formulation of the spectrally dependent flux 
parameterizations. The local fluxes were also calculated in 
the Eddington modification, however, no explicit concatenation 
process occurred. The approach for LOWTRAN 7 is to evaluate the 
local fluxes and then combine them to produce a composite flux 
profile using the adding method. 

In LOWTRAN 7 the appropriate infrared parameterization is 
given by a linear Planck function relationship across an atmos¬ 
pheric layer. This leads to parameterizied two-stream solutions 
for emission from the layer top and layer bottom, and for the 
transmission and reflection functions given explicitly by 


Isaacs et al. (1986): 
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F + - a(PB t -mQ -B fa ) / D (4.42) 

F“ - a(PB, +mQ -B. ) / D (4.43) 

D t 

X - a / D (4.44) 

R - uv(e tl - e zl) / D. (4.45) 


Here the and B, are the full Planckian function for the 
t D 

thermal emission at the layer top and bottom. The individual 
components which make up the expressions above are: 


a - J (1-w ) 

0 

(4.46) 

m - (B - B ) / t 
b t 

(4.47) 

_ rl -tl 

P » ve + ue 

(4.48) 

tl -tl 

! » ve - ue - a 

(4.49) 

„ 2 tl 2 -tl 

(4.50) 

u - (1-a) / 2 

(4.51) 

v - (1+a) / 2 

(4.52) 

tl ■ /3 at. 

(4.53) 


The optical thickness (t) and the single scattering albedo («) 
are parameterized by the following expressions: 

t - ku + t g (l-g) + (4.54) 

« - t: (1-g) / t. (4.55) 

Q S 


Here (k) is the gas absorption coefficient (for a particular 
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wavenumber and probability interval), (u) is the gas density, 
(t ) is the scattering optical thickness, (r ) is the absorp- 

S 3 

tion optical thickness, and (g) is the asymmetry factor for the 
particulate matter in the layer. 

The infrared parameterization given above is significantly 
different from the previous derivation for the thermal fluxes, 
but there are some parallels. The first difference is that my 
method did not use the Planck function explicitly for the top 
and bottom emissions. But, the above method does in fact make 
the same assumption about a linear relationship between the 
Planck funtion and the change in optical depth. Another differ¬ 
ence is that the Eddington modification did not use parameter- 
izied values for the optical depth (thickness) or single 
scattering albedoes-. The Eddington method used the LOWTRAN 
derived transmittances which inherently contain the total 
extinction coefficient based on molecular and particulate 
scattering and absorption extinction coefficients. However, the 
parameterization above does, like my method, use only the layer 
optical depths rather than the total optical thicknesses. It 
is justifiably done in this case because the adding method 
allows layer interaction explicitly while the Eddington modifi¬ 
cation did not relate the adjacent layers. Another difference 
is the use of the exponential terms (e ±r ^) to directly diminish 
or enhance the Planckian emission. My formulation relied upon 
the solution of the matrix vector x to solve for the fluxes at 
each level with the terms containing the source function used 
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to calculate the constant vectors. Thus, no direct scale 
interaction was allowed. Of course, the most significant 
difference lies in the fact that layer transmission and reflec¬ 
tion functions are used to determine the intermediate composite 
fluxes at the layer interfaces. The modification that was 
attempted did not take the reflectance nor transmission value 
into consideration. This did not allow layer interaction which 
caused my flux values to become inflated. No explicit attenua¬ 
tion process was available to diminish the flux values. 

The infrared flux parameterization then represents the 
complex interaction of the effects of molecular and aerosol 
extinction processes related to the thermal emission from an 
isolated' layer. The ratios of single scattering albedo and 
asymmetry factor are used to construct a two-stream approxima¬ 
tion for the amount of longwave radiation being generated from 
layer (N) based on its temperature and optical depth. 


4.3.4 The Adding Method 

The multiple scattering (MS) capability in LOWTRAN comes 
about by combining the stream approximation (SA) with the 
adding method to obtain the required MS flux profiles. Funda¬ 
mentally, the local MS source function is composed of contribu¬ 
tions scattered from layers, both above and below the current 
layer in addition to local thermal emission. Once obtained, the 
source function for the MS radiance may be summed along a 
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desired path to obtain the thermal emission via summation of 
the all local layer Planck function contributions. 

In the adding method, properties of individual layers 
(such as emission, transmission, and reflection) are combined 
algebraically to yield the overall properties of the combined 
layers. A typical application of the adding method consists of 
two opposite unidirectional passes throughout the model atmo¬ 
sphere. The upward pass provides an opportunity to evaluate 
and merge radiance contributions due to thermal emission. In 

addition, the upward pass calculates: (a) layer single scatter- 

N N 

ing albedo (w ), (b) phase function asymmetry factor, (g ), 

which is used to calculate the backscatter fraction, (|5), (c) 

layer reflection and transmission functions R and T t , (d) 

N N 

intrinsic upward and downward fluxes (F*) using the two-stream 

N 

approximation interfaced with the k-distribution method, (e) 
beginning of the adding method for composite upward fluxes 
(F^*) and reflection (R + ), and layer and merged values of the 
thermal radiance contributions based upon the parameterizied 
thermal source function Eqs.(4.42-4.45) 

To illustrate an application of the adding method to cal¬ 
culate the actual flux profiles, consider a layered atmosphere 
shown in Figure (4.3). The upward adding pass begins by adding 
layer (N-l) to layer (N-2). Subsequent layers are added one at 

a time until layer (N+2) is reached. The composite upward flux 

1 ♦ 

from the combination of (N-2) and (N-l), F. t , consists of three 

N 

contributions (see Figure 4.3): 
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Figure 4.3 LOWTRAN 7 Example of the model used for the adding 
method. and T N represent the layer N reflection and 
transmission functions, respectively. F± is the "local" flux 
intrinsic to the layer top and bottom for the infrared parame¬ 
terization. 
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♦ 

i. The intrinsic upward flux from layer (N-l), F N _^ 

ii. The multiple reflection of the intrinsic upward flux from 

♦ 

layer (N-2),transmitted through layer (N-l), F j, 


F N-2 +T N-l + F N-2 T N-1 R N-l R N-2 + F N-2 T N-1 (R N-l R N-2 ) + ' 

F N-2 Vi (1 “ ^-1^-2 - 1 


(4.56) 


where the surface boundary emission is given by: 


F „ - F -ire B(T ) 
N-2 o s s 


(4.57) 


and the reflection funtion for (N-2) is 


R N-2 - R o “ (1 ' C s ) 


(4.58) 


F N-2 ,T N-2 ” F o’ T o “ ° 


(4.59) 


where T is the surface temperature and e is the surface erais- 
s s 

sivity. 


iii. The multiple reflection of the intrinsic downward flux 
from layer (N-l), , reflected from layer (N-2) and 
transmitted through layer (N-l), i.e., 


F N-l T N-l R N-2 + F N-1 T N-l R N-2 (R N-l R N-2^' 
F N-l T N-l R N-2 (1 “ R N-l R N-2 ) ’ 


(4.60) 


These will sum to the composite flux: 


F N-1 * F N-l +T N-l^ F N-2 +F N-l R N-2 )< ' 1 ” R N-l R N-2 ) 


(4.61) 


Similar considerations lead to the composite reflection of 
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Since there are no layers above (N+2), the upward composite 
1 + 

*N+2 1 “ enti ’ call y the actual upward flux, i.e., 


_2* J* 

F * F 
N*2 N*2 


(4.66) 


As is evident from the above discussion, calculated compo¬ 
site layer values are required for the evaluation of composite 
values of subsequent layers. This process was not used in the 
Eddington modification. Omission of the layer connections 
became the fatal physical flaw since no accountability of layer 
interaction was included. 
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At the termination of the upward pass, the following quan¬ 
tities would have been calculated for each layer and stored for 

N 

the downward pass: (a) the intrinsic layer values, u , R , T, 

on N 

— 1 * + 

F v _ , (b) the composite quantities, F VT , R , and the non- 

N N 

scattered contribution to emission, E^, where 

N 

E N - E N ' (1 'V J SA (*•«> 



(4.70) 






56 


At this point, actual upward and downward fluxes are available 
at the boundaries of layer (N+2). This approach to secure the 
actual upward and downward fluxes which is radically different 
from the Wiscombe (1977) procedure where flux continuity is 
assumed at the layer interfaces i.e., the components and H. 
contributing to the upward and downward fluxes are equal. No 
distinction is made concerning the optical properties of the 
adjacent layers, whereas, LOWTRAN 7 does distinguish those pro¬ 
perties based on the reflectance, transmittance, and composite 
fluxes. 

2 * 2 - 

The actual (F , upward and downward fluxes are 

used to calculate the multiple scattering source function 

MS 

applicable to that level e.g., us * n § t * ie stream approxi¬ 

mation Eq. (4.39).. This source function provides the layer 

MS 

contribution to the multiply scattered emission, E. t „, and this 
is added to the thermal contribution, Eq- (4.67) 

evaluated during the upward pass. Proceeding downward, the com¬ 
posite quantities are first calculated by adding the next layer 
down to the composite above it, for example, the downward com¬ 
posite flux emerging from layer (N+l) and the corresponding 
downward composite reflections are: 

F N*l* F N*l +T N+l (F N-2 + F N*1 R N*2 )(1 " ^+2 

(4.71) 

^♦1 " ^+1 + ~ ^+1 ^* 2 ' ' 

Actual fluxes at the interface of layers (N) and (N+l) are then 
obtained from the composite values: 
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4 ‘ ■ ( 4 ' * 4 m 4 )(i - 4 4 - i >' 1 

2- I— i ♦ — ♦ — — I 

F„ - ( F„ . + F„ R f ,)(1 -R R ,) . 

N*1 N*1 N N*l /V n N*1 


(4.72) 


For the next layer, N, the procedure is identical, resulting in 
the composites: 


4‘ ■ 4 *v4m * 44.i )(1 - 4.i rl 

4 - ■'h * Vi 4 (1 - “n Vi>' 1 - 


(4.73) 


and actual fluxes: 


4-i ■ < 44 * 4' 4-i )(1 - 4-i 4 r ‘ 
4' - <4'* 44 4 )(i - 4-i 4 )_1 - 


(4.74) 


Finally, the procedure terminates at the surface layer (N-2), 
(N-l) interface where: 


44 ■ 4-i * t h-i< 4~ * 4-i4 )(1 - Vi 4>"‘ 


b n-i " r k-l + Vh-i " 1 " r n-: r n j 


(4.75) 


Then the total fluxes become: 


F N -2 ” ( F N -2 + F N -1 R N- 2^ 1 ~ V 2 VP 

F N-i “ (F i-l + F N -2 Vl )(1 - *N -2 VP' 1 * 


(4.76) 


♦ 

In this last step, the upward surface emission flux, F^ ^ and 
reflection function, , are composite quantities since there 
are no layers below. This completes the downward pass and the 
multiply-scattered radiance calculation. 
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4.3.5 Path Radiance 


Path radiance calculations are generated with the combina¬ 
tion of the stream approximation (SA), and the multiple 
scattering source funtion at optical depth, (c), for a path 
zenith angle cos For level N, the level source function is 
given by: 


N(±p) 


u-“ 0 ib<v> 


(4.77) 


N , r _±l - .N. n. ± N , N. ■, 

+ « / n [F Mpg ) + F 8 (H g )]• 

O N N 


The first term is due to the non-scattered thermal emission, 
while the second terra is due to multiple scattering. Thus the 
source function can be written as 


t SA T o . N _ n t MS 
J N = (U o ,T N ) + J N 


(P, v g N . 


(4.78) 


where the thermal emission term J° depends only on local layer 

MS 

properties, and the multiply-scattered term requires, in 

addition to local properties, the local fluxes which depend on 
the overall path properties. 


The path radiance or emission with MS is obtained by sum¬ 
ming the radiance contributions along a particular path. 

Evaluation of the multiply-scattered radiance contribution, 
MS 

E._ , for each layer is accomplished via the adding method, as 
N 

explained above, on the complementary downward pass where it is 
merged with the previously calculated non-scattered contribu¬ 
tions to obtain the total radiance. Once obtained the radiances 







59 


I + and I will be multiplied by n to produce the isotropic 
hemispheric flux F* and F+ for the respective level. 

4.3.6 K-distribution Method 

The basic problem encountered in the calculation of radia¬ 
tive transfer in low spectral resolution models like LOWTRAN, 
in the prescence of hazy or cloudy conditions, is the coupling 
between the processes of scattering and absorption by atmos¬ 
pheric gases and that of cloud/aerosol particles. It is 
extremely difficult to distinguish the two attenuation effects 
from one another. The complicated molecular absorption band 
structure for gases would not allow a rapid line by line fre¬ 
quency integration. However, making certain assumptions about 
the line shape one can calculate from available data the 

absorption coefficient for gases, (k (v)), versus wavenumber 

a 

for given temperatures and pressures. 

In principle, one can directly employ k^Cv) in the radia¬ 
tive transfer equations and obtain the numerical solutions. As 
pointed out by Arking and Grossman (1972), the only difficulty 
is the complicated nature of the function (k (v)). A typical 

cl 

band will contain anywhere from several hundred to several 
thousand lines that must be included to obtain reasonable accu¬ 
racy . 

A solution to this difficulty can be found through the use 


of the k-distribution function. This method makes use of the 
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fact that for a homogeneous atmospheric layer, the transmission 
within a relatively wide spectral interval is independent of 
the ordering of the value of k^ with respect to wavenumber, but 
depends upon the fraction of the interval that is associated 
with a particular value of k^. 

In a homogenous gas layer, the k-distribution function is 

formally related to the mean transmission function T (u), 

Av 

_ , „ , ,. r -ku dv 

T. (u)= 1/Av J e 
Av 4 

Av 

- J f(k)e' kU dk 
o 


-ku 

e dg 


(4.79) 


o 


n -k.u 
- Ee 1 Ag. 

1 

where (Av) is the narrow repeated interval (20 cm * in LOWTRAN) 
and (u) is the gas amount. The f(k) for a given gas at a speci¬ 
fied Av is the probability density function such that f(k)dk is 
the fraction of the frequency interval for which the absorption 
coefficient is between k and k+dk. Eq. (4.80) reveals that the 
transmission depends on the distribution of the k-values within 
Av, but not on the ordering of the values. 


The cumulative k-distribution function is g(k) while 
(k^,Ag^) are the discrete sets of values to approximate the 
integral in Eq. (4.80). For lines which are randomly distri¬ 
buted, allowing overlap, and assuming the lines shapes are 
square, the k-distribution function can be a Poisson distribu¬ 
tion according to Arking and Gossraan (1972). For example, if 
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one line contributes an absorption k Q over the interval 2ad 
where (d) is the mean line spacing, then the permitted values 
are k n k Q , where n is a positive integer, and the k- 
distribution function is 

g(k)(k=nk o ) - 2(a) n e' 2a/n! , n-0,1,2. (4.80) 

Transmission functions for most atmospheric gases and 
perhaps also at almost all spectral ranges are non-exponential, 
if derived for spectral intervals of finite width AX.. The 
theoretical treatment of radiative transfer in scattering and 
absorbing media, however, requires strict exponential transmis¬ 
sion laws. 

This deficiency is overcome by fits of the transmission 

function (T) to finite sums of exponentials, as defined by 

_ n -b.u 

T a (v) = Z a. e 1 , (4.81) 

Av . , l 

1 = 1 

where v is the center wavelength of the interval Av, and u is 
the absorbing path length. To obtain physically meaningful 
terras of this sum, the coefficients a^ and exponents b^ must 
obey the following conditions: 

n 

0 < a. < 1 , Z a. = 1 
1 i-1 1 

and 

b. > 0. 
i 

The exponents b^ are essentially identical with the volume 
absorption coefficients (k) and have the dimension of (km *). 

By expressing the band model transmission as the sum of 
exponentials, the multiple scattering calculation for each 
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component can be performed independently as if it were a mono¬ 
chromatic problem. Recall, that the optical depth in the 
infrared parameterization was related to the extinction coeffi¬ 
cient by: 

t - k u + z (1—g) + t . (4.82) 
v s a 

This effectively allows differentiation between the aerosol 
contribution to extinction and gaseous effects of absorption 
and scattering. 

The k-distribution method decouples the band model deter¬ 
mined optical properties into a set of equivalent monochromatic 
calculations which are summed to give the spectrally averaged 
results. For LOWTRAN 7, the stream approximation is performed 
through an interface routine consisting of the k-distribution 
method. This combination of methods, namely the stream approxi¬ 
mation and the k-distribution method have been used to study 
the radiative effects of aerosols (see Hansen et al., 1980). 
My Eddington modification did not refine the difference between 
the individual absorption bands for gases but took the extinc¬ 
tion coefficients from the inverted transmittances from LOWTRAN 
6 . 

The fitting routine of Wiscorabe and Evans (1977) has been 
used for the two LOWTRAN transmission functions of water vapor 
/uniformly mixed gases and ozone the primary gaseous absorbers 
in the troposphere. The accuracy of the fitting is, in general, 
within a few percent for transmittances greater than 0.1. For 
inhomogeneous atmospheres, the adopted scaling approximation is 
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where T , 
o 

tively. 


k.(P.T) - k i (P 0 ,T o ) P/B o AT 0 /T) 

P are reference temperature and pressures, 
o 


(4.84) 

respec- 





5.0 PHYSICAL PARAMETERS for the AEROSOL-LADEN 
BOUNDARY LAYER 

5.1 Overview 

Armed now with a code which can adequately handle the com¬ 
plexity of multiple scattering, the comparison between multiple 
and single scattering can proceed. The method of inspection is 
to use the principle of juxtaposition to elicit the quantified 
difference between the net fluxes at prescribed levels in a 
heavily aerosol-laden boundary layer representative of a natur¬ 
ally occurring haze condition. In my experience, this occurs 
during the summer months under high pressure conditions with a 
well established inversion trapping aerosols within the lower 
troposphere. The maximum input of aerosols may be achieved in 
the vicinity of urban areas where anthropogenic materials are 
added to the background distribution. 

5.2 Model Assumptions 

To accomplish my topical objective of quantifying the 
difference between multiple and single scattering in an 
aerosol-laden boundary layer, a model representative of an 
urban, summer time scenario is selected. Such a model is avail¬ 
able in the LOWTRAN family of atmospheric models The overall 
structure is that of a plane-parallel, cloudless environment 
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supplemented with a tropical water vapor profile. This model 
atmosphere is subdivided into 11 levels and ten layers starting 
at two kilometers and extending downward to the surface which 
is assumed at sea level. Table 5.1 shows the profile of various 
parameters applicable to the following discussion. 

Several simplifying assumptions are made concerning the 
model and the method of comparison. 1) Scattering is done only 
by aerosols. Molecular scattering is ignored. The justification 
for this lies in the fact that gaseous scattering (Rayleigh 
scattering) is inversely proportional to the fourth power of 
the wavelength. Therefore, in the terrestrial infrared domain 
(£ 5.0pm) the effect is negligible. 2) Particle temperature is 
equal to the medium temperature. 3) Kirchhoff's law applies. 
This law states a medium may absorb radiation of a particular 
wavelength, and at the same time also may emit radiation of the 
same wavelength. The rate at which emission takes place being a 
function of temperature and wavelength, under the condition of 
local thermodynamic equilibrium (LTE). In other words, the 
reduction in the intensity of a pencil of radiation is a conse¬ 
quence of absorption and scattering by the intervening 
material. However, where absorption is occurring so is emission 
according to Planck's law of blackbodies. Therefore, the inten¬ 
sity may be strengthened by emission within the line of sight 
plus multiple scattering from all directions into the line of 
sight. Symbolically, the above physical processes are 


represented by: 







TABLE 5.1 BOUNDARY LAYER MODEL ATMOSPHERE 


L 

Z 

RH 

P 

T 

VIS 

V 


(km) 

(Z) 

(mb) 

(°K) 

(km) 

(km ) 

1 

2.000 

100.0 

802.0 

289.7 

0.55 

7.06 

2 

1.750 

97.3 

826.2 

286.4 

0.83 

4.78 

3 

1.S00 

95.2 

850.9 

287.6 

1.15 

3.41 

4 

1.250 

92.9 

876.1 

288.6 

1.53 

2.51 

5 

1.000 

91.1 

920.0 

289.6 

1.94 

2.01 

6 

0.750 

88.9 

928.6 

290.8 

2.46 

1.59 

7 

0.500 

88.2 

955.9 

291.9 

2.96 

1.32 

8 

0.250 

87.3 

984.1 

293.0 

3.47 

1.12 

9 

0.100 

86.7 

1001.3 

293.7 

3.79 

1.01 

10 

0.050 

86.4 

1007.1 

294.0 

3.89 

0.99 

11 

0.001 

86.2 

1012.9 

294.2 

3.99 

0.97 


Z - elevation of level L in (km) 

RH - relative humidity at level L 

P - pressure at level L 

T - temperature at level L 

VIS - meteorolgical range (3.912/E t ) 

E t » extinction coefficient norraalizied at 0.55pm 
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j - k B(T) (5.1) 

where j is the source function coefficient and k is the mass 
v v 

extinction cross section (the sum of the mass absorption and 
scattering cross section) both coefficients having the same 
physical meaning. 4) The visibility in my model decreases with 
height to reach its minimum at the top of the inversion layer. 
Correspondingly, the extinction coefficient increases to reach 
its maximum at the same location. This model adjustment is 
selected to mimic summertime conditions in an urban environment 
where, in the presence of an inversion and strong insolation, 
convective "boiling" and vertical turbulent transport differen¬ 
tially mix the particulate matter with lowest visibilities 
found immediately below the inversion, as described in a model 

by Venkatram and Viskanta (1977). 5) The number concentration 

4 3 

is fixed at 2.0x10 particles/cm . However, particle size 
increases with increasing humidification. This is a LOWTRAN 
artifact which relates the meteorological range to the layer 
relative humidity based on extensive research conducted by 
Hanel (1972, 1976) and McClatchey et al. (1972). dependent 

configuration. It is intented only to quantify the difference 
of multiple and single scattering allowing only the transfer of 
radiant energy to occur over the same environmental conditions. 

5.3 Urban Aerosol Characteristics 

Urban atmospheric aerosols and their radiation absorption 
and scattering properties are characterized by their size dis- 
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tribution, chemical composition, and physical properties such 
as density and the complex index of refraction. The urban aero¬ 
sol can be separated into distinct groups based on their size 
and chemical composition according to Viskanta and Weirich 
(1979). The first group consists generally of particles with 
radii greater than 2pm. These particles are composed primarily 
of silicon oxides with substantial amounts of calcium and 
lesser amounts of iron oxides and metal chlorides. The second 
group consists generally of particles which have radii less 
than 2pm. These particles are chiefly anthropogenic in origin 
and are either produced directly as aerosol or by heterogeneous 
nucleation of gas-phase pollutants. Chemical analysis of these 
particles shows large concentrations of sulfur and carbon com¬ 
pounds with varying amounts of lead, zinc, arsenic, and other 
metals. 

The urban aerosol model in LOWTRAN 7 is a combination of 
the rural aerosol, composed of a mixture of 80 percent of water 
soluble substance (ammonium and calcium sulfate and also 
organic compounds), and 20 percent sootlike, carbonaceous aero¬ 
sols. The sootlike aerosols are assumed to have the same size 
distribution as the other components according to Shettle and 
Fenn (1979). Because of the hygroscopic nature of the many 
chemical compounds which make up the urban aerosol and the 
presence of high humidities in urban areas at mid-latitudes in 
the summer, the particles are often composed largely of liquid 
water. This then leads to water coated particles which are com- 
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posite spheres containing several individual constituents. Fig¬ 
ure (5.1) shows a visual presentation of this process. Since 
the volume increases by the cube of the radius, the fraction of 
water also increases rapidly when the air humidity increases. 
This means that the Mie theory for spherical particles will be 
more accurate the higher the air humidity is and also that the 
refractive index (real portion) approaches the value of liquid 
water (1.33). 

5.4 Humidification Effects on Aerosol Properties 

As indicated above, the urban aerosols being hygroscopic 
in nature are subject to changes in their physical properties 
as the ambient relative humidity increases. Water vapor con¬ 
densing out of the atmosphere onto particulates suspended in 
such an atmosphere tends to swell their size and change their 
composition thus their effective complex refractive index. The 
resulting effect of the aerosols on the absorption and scatter¬ 
ing of light will correspondingly be modified. 

The change in particle size is related to the relative 

humidity by an equation developed by Hanel (1976): 

r(a ) ■ r [1 + p ra (a /m .. 

w o w w o (5.2) 

where r is the dry particle radius, p is the particle density 
o 

relative to that of water, ra is the dry particle mass, m (a ) 

’ o w w 

is the mass of the condensed water, and (a ) is the water 

w 

activity which is essentially the relative humidity (f), 








1. Dry porticlt (RH«0%) 



water soluble substonces, density s 6 S 
water insoluble substances, density = J, 
radius * r a 

fraction ol water soluble substances 


2 Droplet (RH» 8S%) 


-salt solution, density »6 

--woter soluble substances, denstty*5 s 

--woter insoluble substances, density * 6 , 

-radius = r 

dissolved fraction of E - C 


3 Droplet (RH>95%) 



Figure 5.1 An aerosol particle model which shows an increase 
in radius (r) and change in the fraction of water soluble and 
water insoluble substances (e) as the relative humidity 
increases in the environment refelcted in the densities 5. 
(after Nilsson, 1979). 1,8 
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corrected for the curvature of the particle surface, vis., 

a w - f exp(-2oV w /R w T r), (5.3) 

where o is the surface tension on the wet particle surface, V 

w 

is the specific volume of water, R^ is the specific gas con¬ 
stant for water, and T is the absolute temperature. 


Once the wet aerosol particle size is found from Eqs (5.1 

and 5.2), the effective complex refractive index for composite 

water coated substances becomes the volume weighted average of 

the refractive indexes of the dry aerosol substance, n and 

o 

water n^. Equivalently, this can be written as (after Hanel, 


1976): 


n = n 


w 


o w o w 


(5.4) 


Figure (5.2) shows the real and imaginary parts of the complex 
refractive index as a function of wavelength for an urban aero¬ 
sol collected in a study by Nilsson, (i979). The real part 
decreases when the humidity increases except around 3.15pra and 
8.2pra. This is particularly noticeable for wavelengths greater 
than 9.0pra. The imaginary part has a pronounced peak at 9.2pm 
caused by absorption in quartz, silicon, and ammonium sulfate. 
This peak decreases when the air humidity increases simultane¬ 
ously as peaks caused by absorption in water, grow at 2.95pm, 
6.0pm, and 15-20pm. Greater absorption of electromagnetic 
energy by the water substance (vapor and liquid) in the above 
opaque spectral bands is indicated by the increasing values of 
the imaginary portion of the complex index of refraction. Thus 
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Figure 5.2 The real and imaginary parts of refractive index 
for aerosol partilces as a function of wavelength at different 
relative humidities for an urban aerosol collection (after 
Nilsson, 1979). 
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more energy will internalized by the intervening water vapor 
and water coated aerosols. This reduces the probability that 
the scattering and refraction processes will contribute to the 
diffuse radiation field. However, increased absorption leads to 
a concurrent greater emission from the intervening medium. 
Therefore, the "true creation" of photons, a direct contribu¬ 
tion, would enhance the total flux density field. These physi¬ 
cal phenomenon are intertwined in a complex manner in a model 
where the visibility is tied to the relative humidity. The 
change of the complex refractive index along with the increas¬ 
ing size of the particles as the humidity approaches saturation 
effectively changes the extinction coefficient (the combination 
of scattering and absorption bulk properties described in 
Chapter 3). This in turn effects the meteorological range (V) 
given by the Koschmieder formula 

V - 1/0 In 1/e - 3.912/0 (5.5) 

where (0) is the extinction coefficient, and (e) is the thres¬ 
hold contrast, set equal to 0.02. As used in LOWTRAN, the 
inputs for visibility are in terras of meteorological range, 
with (0) evaluated at 0.55pm. 

The meteorological ranges presented in Table (5.1) come 
from the Army Vertical Structure Algorithm (VSA) build into 
LOWTRAN 7. This is a subprogram which effectively subdivides 
the model boundary layer into a finer resolution. The VSA 
allows various prescriptions for the aerosol configurations to 
include the heavy haze condition (case 2) which forms the basis 
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of ray heavily aerosol-laden test case. 

5.5 Relative Fractional Comparisons and Net Flux Generation 

LOWTRAN 7 is run in the radiance mode to calculate the 

difference between multiple and single scattering. This is 

accomplished by toggling between the two available modes 

currently installed. The method of comparison is then to find 

the net fluxes (F ) for each scattering configuration at set 
net 

levels and compare the relative fractional difference (reldiff) 
given by: 

(MS-SS)/SS ( 5 - 6 ) 

where MS is the multiply scattered net flux values and SS is 
the singly scattered net flux. 

To calculate the net fluxes (F ) at each level the 

net 

respective upward and downward radiances (If and I*) must be 
known. The net flux ( F net > is related to these by the simple 
relationship 

F - ir (It - 1 + ) (5.7) 

net 

based on the isotropy conditions in the thermal radiation 
fields. Here the downward radiance (I*) is chosen as subtrahend 
since it is expected that the upward radiance at each level 
will be larger since the highest temperatures are located in 
the lowest levels of .. v boundary layer. Thus, emission would 
be larger looking downw, 
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The upward radiances are calculated by looking down toward 
the surface from each successive preestablished level starting 
at two kilometers. The surface boundary conditions are a sur¬ 
face temperature(T ) of 300°K and surface emissivity (c) of 

s 

0.85. These values are selected to mimic mid-latitude summer¬ 
time conditions and a combination of rural and urban landscapes 
as suggested by Goody (1964). 

The LOWTRAN 7 program is set to move through the levels 
once for each of 34 wavenumber intervals beginning at 350 cm ^ 
to 2000 cm * at 50 cm * steps. This is done to generate a spec¬ 
trally sensitive output which reflects the absuiption and emis¬ 
sion characteristics of the radiatively active constituents in 
the test model. An angular dependence is set by looking down at 
zenith angles of 180°, 140°, and 120°. 

The downward radiances (I*) are calculated by looking up 
toward the top of the boundary layer (2 kms) with the exception 
that the first downward radiance is a composite of the inter¬ 
vening layers from 20.0 kms. Successive downward radiance cal¬ 
culations occur at each level in the boundary layer again vary¬ 
ing the zenith angle, looking up, at 0°, 40°, and 60° then over 
same spectral intervals. 

Each respective level up and down radiance is then multi¬ 
plied by n to produce the hemispheric isotropic flux. Assuming 
isotropy for these fluxes is consistent with the type of radi¬ 
ant emission from layers under LIE. This is again a basic 
assumption pertinent to my study. Now that the fluxes are 
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available Eq. (5.5) can be used to calculate the net fluxes 


(F 


net 


)• 


The relative fractional difference is produced using the 
net fluxes at each level for each type of scattering mode (SS 
or MS) with comparisons being made over zenith angle (9) and 
wavenumber interval (Av). The formulation appears as: 

we e e 

Reldif(9,Av) » [Fnet (9,Av) - Fnet (0,Av)J 

ss (5 ' 8) 

/ Fnet (9,Av). 


The information from Eq. (5.7) provides insight about the quan¬ 
tified difference between MS and SS in a heavily aerosol-laden 
boundary layer environment. 




6.0 Interpretation of Results 


6.i Overview 

Careful examination of the output data revealed several 
distinct trends in regards to the behavior of the relative 
fractional difference (Reldiff) and particular spectral 
regions. It is evident that the distinction between multiple 
scattering (MS) and single scattering (SS) is practically zero 
for those wavenumber regions where strong water vapor absorp¬ 
tion is occurring. Uniformly, over all zenith angles, this 
trend repeated itself at elevations in the boundary layer where 
the relative humidities exceeded 88Z (level 7, 0.50 km). In 
these same spectral regions it is found that, below level 7, 
negative MS net fluxes (Ft - Ft), downward directed net flux, 
dominate by as much as 4-0 times the SS net fluxes. Another 
finding is that for "window regions" in the infrared portion of 
the spectrum, which are much less opaque to thermal radiation, 
the MS net fluxes exceed the SS net fluxes at all levels and 
increases to a maximum as the zenith angle increases. There¬ 
fore, since the above trends are found throughout the data, the 
presentation for the interpretation of results will be limited 
to two representative spectral bands, as numerically exempli¬ 
fied by the 1650 cm * (6.1pm) and 950 cm * (10.5pm) 
wavenumbers. 


77 







78 


6.2 Atmospheric Effect and the Thermal 
Infrared Spectrum 

Certain portions of the infrared radiation spectrum are 
absorbed by various gases in the atmosphere. Among these gases, 
carbon dioxide, water vapor, and ozone are the most important. 
Figure 6.1 shows the emission spectrum for a surface tempera¬ 
ture of 295°K. Strong absorption is found due to water vapor in 
the 6.3pm band from about 1200 to 2000 cm * and in the rota¬ 
tional bands below 600 cm Except for ozone, which has an 
absorption band in the 1040 cm * (9.6pm) region, the atmosphere 
is relatively transparent from 800 to 1200 cm *. This region is 
referred to as the "atmospheric window." It is seen from Figure 
6.1 that the choice for the 1650 cm 1 band places it squarely 
in the middle of the strongest H^O absorption. Whereas, the 
950 cm ^ band is selected to be outside the ozone (0^), carbon 
dioxide (CO , and water vapor (H^O) bands. 

6.3 The 1650 cm ^ (6,0pm) Absorption Band 

The 1650 cm ^ absorption band comparisons of MS versus SS 
fluxes are shown in Table 6.1 through 6.3. These results show 
that regardless of the viewing angle above 0.25 km (levels 1 
through 7) the difference between MS and SS net fluxes are 
eradicated. My interpretation of this result is that water 
vapor plus liquid water coated aerosol absorption is so great 
that any distinction between the scattering processes are com¬ 


pletely negated. 





Figure 6.1 The terrestrial infrared spectra and various 
absorption bands. Also shown is an actual atmospheric emission 
spectrum taken by the Nimbus IV IRIS instrument near Guam at 
15.1 N and 215.3 W (after Liou, 1980). 
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TABLE 6.1: 1650 cm" 1 ABSORPTION BAND 
ZENITH ANGLES 0 # and 180“ 


L 

SS 

MS 

Reldiff 

1 

0.9400E-07 

0.9400E-07 

0.0000 

2 

0.3100E-07 

0.3100E-07 

0.0000 

3 

0.6300E-07 

0.6300E-07 

0.0000 

4 

0.9400E-07 

0.9400E-07 

0.0000 

5 

0.9400E-07 

0.9400E-07 

0.0000 

6 

0.9400E-07 

0.9400E-07 

0.0000 

7 

0.1260E-06 

0.1260E-06 

0.0000 

8 

0.1570E-06 

-.5027E-05 

-33.02 

9 

0.1260E-06 

-.5089E-05 

-41.39 

10 

0.1570E-06 

-. 5088E-05 

-33.41 

11 

0.1250E-06 

-.2526E-05 

-21.21 


L : level number 
SS 

SS *■ Fnet : the single scattered net flux at level L 
MS 

MS - Fnet : the multiple scattered net flux at level L 
Reldiff : relative fractional difference between the 
net fluxes (MS-SS) / SS 
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TABLE 6.2: 1650 cm" 1 ABSORPTION BAND 
ZENITH ANGLES 40° and 140° 


L 

SS 

MS 

Reldiff 

1 

0.9400E-07 

0.9400E-07 

0.0000 

2 

0.310QE-07 

0.3100E-07 

0.0000 

3 

0.6300E-07 

0.6300E-07 

0.0000 

4 

0.9400E-07 

0.9400E-07 

0.0000 

5 

0.9400E-0 7 

0.9400E-07 

0.0000 

6 

0.9400E-07 

0.9400E-07 

0.0000 

7 

0.1260E-06 

0.1260E-06 

0.0000 

8 

0.1570E-06 

-. 5027E-05 

-33.02 

9 

0.1260E-06 

-.5089E-05 

-41.39 

10 

0.1570E-06 

. 5089E-05 

-33.41 

11 

0.1250E-06 

1885E-05 

-16.08 


L : level number 
SS 

SS = Fnet : the single scattered net flux at level L 
MS 

MS =* Fnet : the multiple scattered net flux at level L 
Reldiff : relative fractional difference between the 
net fluxes (MS-SS) / SS 
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TABLE 6.3: 1650 cm" 1 ABSORPTION BAND 
ZENITH ANGLES 60* and 120 # 


L 

SS 

MS 

Reldiff 

1 

0.6300E-07 

0.6300E-07 

0.0000 

2 

0.3100E-07 

0.3100E-07 

0.0000 

3 

0.6300E-07 

0.6300E-07 

0.0000 

4 

0.9400E-07 

0.9400E-07 

0.0000 

5 

0.9400E-07 

0.9400E-07 

0.0000 

6 

0.9400E-07 

0.9400E-07 

0.0000 

7 

0.1260E-06 

0.1260E-06 

0.0000 

8 

0.1570E-06 

-. 5027E-05 

-33.02 

9 

0.1260E-06 

-. 5089E-05 

-41.39 

10 

0.1570E-06 

-. 5089E-05 

-33.41 

11 

0.1250E-06 

1884E-05 

-15.83 


L : level number 
„ SS 

SS = Fnet : the single scattered net flux at level L 
MS 

MS « Fnet : the multiple scattered net flux at level L 
Reldiff : relative fractional difference between the 
net fluxes (MS-SS) / SS 
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Sensibly, this incorporates the notion that as the rela¬ 
tive humidities increase the size of the aerosol particles 
would increase. Therfore, an increasing amount of liquid water 
is needed to encapsulate the originally drier particle. This in 
turn changes the complex index of refraction to higher ima¬ 
ginary values as pointed out in Figure (5.2). Appropriately, 
this allows more absorption to occur within those layers where 
the relative humidities are high enough to trigger the effect. 
This is graphically depicted in Figure 6.2 which is a plot of 
the partial aerosol extinction coefficients (absorption and 
scattering) at 50% and 99.8% relative humidities (RH) for an 
urban aerosol which is similar in number density, composition, 
and size distribution used in my investigation. At the 99.8% 
RH at wavenumber 6.0pm the scattering and absorption coeffi¬ 
cients are approximately equal. However, closer inspection 
shows that the scattering extinction coefficient (Eg) for this 
spectral band is decreased and the absorption extinction coef¬ 
ficient (E ) increases to a maximum around 15.0pm. This 
a 

shows, indirectly, that the probability increases that a photon 
emitted or scattered (by either process) would be almost 
instantaneously absorbed. In addition, analysis of the up and 
down fluxes for each angle showed that the differences between 
them is minimal for both MS and SS (10 7 ) at levels where the 
RHs exceeded 88% (see Table 6.4 for a representative example). 
This indicates that isotropic thermal emission (evenly distri¬ 
buted up and down fluxes) is occurring in those layers and that 
any directionality due to scattering, implicitly present in the 













TABLE 6.4: 1650cm" 1 REPRESENTATIVE MULTIPLE 
and SINGLE SCATTERING UP and DOWN FLUXES at 
ZENITH ANGLE 40° 


Multiple Scattering 


Level 

down-flux 

up-flux 

net-flux 

1 

0.4304E-05 

0.4398E-05 

0.9400E-07 

2 

0.4210E-05 

0.4241E-Q5 

0.310QE-07 

3 

0.4335E-05 

0.4398E-05 

0.6300E-07 

4 

0.4461E-05 

0.4555E-05 

0.9400E-07 

5 

0.4618E-05 

0.4712E-05 

0.9400E-07 

6 

0.4775E-05 

0.4869E-05 

0.9400E-07 

7 

0.4838E-05 

0.4964E-05 

0.1260E-06 

8 

0.5027E-05 

0.2099E-18 

-.5027E-05 

9 

0.5089E-05 

0.3770E-12 

-.5089E-05 

10 

0.5089E-05 

0.1784E-09 

-.5089E-05 

■ 11 

0.5121E-05 

0.3236E-05 

1885E-05 


Single Scattering 


Level 

down-flux 

up-flux 

net-flux 

1 

0.4304E-05 

0.4398E-05 

0.9400E-07 

2 

0.4210E-05 

0.4241E-05 

0.3100E-07 

3 

0.4335E-05 

0.4398E-05 

0.6300E-07 

4 

0.446 IE-05 

0.4555E-05 

0.9400E-07 

5 

0.4618E-05 

0.4712E-05 

0.9400E-0 7 

6 

0.4775E-05 

0.4869E-05 

0.9400E-07 

7 

0.4838E-05 

0.4964E-05 

0.1260E-06 

8 

0.5027E-05 

0.5184E-05 

0.1570E-06 

9 

0.5089E-05 

0.5215E-05 

0.1260E-06 

10 

0.5089E-05 

0.5246E-05 

0.1570E-06 

11 

0.5121E-05 

0.5246E-05 

0.1250E-06 
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backscatter fraction (0) for MS (0.15 from Figure 4.2) and the 
asymmetry factor (g) (.8520 from Shettle and Fenn, 1979), is 
extinguished. The overall effect is rapid re-absorption and 
re-emission which masks out any contribution by aerosol 
scattering. 

The sudden appearance of the MS net downward fluxes at 

levels eight through eleven (present in Tables 6.1-6.3) is 

explained in conjunction with the ideas presented above. In 

these cases, the MS scattering process supercedes absorption 

and establishes a directional character. So much so as to 

overwhelm the SS process by 40 times! This notion is supported 

by Figure 6.2 with emphasis now shifted to the 50Z RH graph. In 

this case the scattering extinction coefficient (E ) is much 

s 

larger for 6.0ym than the absorption coefficient ( E a )- This 
indicates a trend whereby as the RHs decrease, even marginally, 
the MS process with its favored downward directionality is able 
to produce much larger flux densities due to the presence of 
aerosols. The fact that those layers above, when combined 
together, are by nature much thicker ensures a much larger out¬ 
put of photons in the downward direction. This idea is sup¬ 
ported again by inspection of the downward and upward fluxes 
for each process. In each angular case the MS down fluxes 
overwhelmed the up fluxes by eight orders of magnitude at all 
levels greater than seven (see Table 6.4). Here it can be seen 
that the difference of the SS fluxes do not change appreciably 
showing conclusively the directional nature the MS flux. 
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6.4 The 950 cm ^ (10.5pm) Window Region 

Tables 6.5 through 6.7 present the net flux differences 
for the 950cm ^ window region. In this spectral domain, absorp¬ 
tion by the principal infrared absorbing gases is greatly 
reduced. Therefore, a comparison of MS and SS fluxes can be 
presented with absorption considerations minimized. Several 
trends are noted: First, positive upward net fluxes dominated 
over negative downward net fluxes for all zenith angles and 
elevations. Second, multiple scattering (MS) produces larger 
positive net fluxes than single scattering (SS) over all eleva¬ 
tions with an increasing difference between the two being found 
as the zenith angle (9) increases. Third, for both SS and MS, 
the positive flux densities increase as the level elevation 
from the surface becomes greater. Fourth, a downward trend is 
seen for positive flux densities for both MS and SS as the 
geometric path length expands due to secant elongation. 

The positive upward fluxes are expected, now, since in 
this spectral regime thermal emission for terrestrial tempera- 

r 

tures is approaching a maximum and such emission from lower 
elevations at higher temperatures is capable of being transmit¬ 
ted to the levels above. A cumulative effect occurs as the 
concatenation of levels takes place. Thus, a straight downward 
looking sensor would "see" larger flux density values as the 
distance between the surface and the sensor increased. This is 
similar to what happened with the immense downward flux found 
in the strong absorption band. The preponderance of emitted 




TABLE 6.5: 950 cm" 1 WINDOW BAND 
ZENITH ANGLES 0“ and 180* 


L 

SS_ „ 

Fnet 

MS_ 

Fnet 

Reldiff 

1 

0.2532E-04 

0.2670E-04 

0.0545 

2 

0.2413E-04 

0.2695E-04 

0.1169 

3 

0.2143E-04 

0.2485E-04 

0.1596 

4 

0.1935E-04 

0.2321E-04 

0.1995 

5 

0.1766E-04 

0.2177E-04 

0.2327 

6 

0.1618E-04 

0.2030E-04 

0.2546 

7 

0.1480E-04 

0.1923E-04 

0.2993 

8 

0.1341E-04 

0.1505E-04 

0.1223 

9 

0.1260E-04 

0.1596E-04 

0.2667 

10 

0.1232E-04 

0.1643E-04 

0.3263 

11 

0.1200E-04 

0.1687E-04 

0.4058 


L : level number 
SS 

SS » Fnet : the single scattered net flux at level L 
MS 

MS ■ Fnet : the multiple scattered net flux at level L 
Reldiff : relative fractional difference between the 
net fluxes (MS-SS) / SS 
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TABLE 6.6: 950 era" 1 WINDOW BAND 
ZENITH ANGLES 40° and 140“ 


L 

SS 

MS 

Reldiff 

1 

0.2441E-04 

0.2554E-04 

0.0463 

2 

0.2287E-04 

0.2567E-04 

0.1224 

3 

0.1964E-04 

0.2313E-04 

0.1777 

4 

0.1725E-04 

0.2117E-04 

0.2272 

5 

0.1536E-04 

0.1954E-04 

0.2721 

6 

0.1373E-04 

0.1807E-04 

0.3161 

7 

0.1219E-04 

0.1684E-04 

0.3815 

8 

0.1075E-04 

0.1172E-04 

0.0902 

9 

0.9860E-05 

0.1300E-04 

0.3185 

10 

0.9520E-05 

0.1357E-04 

0.4254 

11 

0.9230E-05 

0.1429E-04 

0.5482 


L : level number 
SS 

S3 - Fnet : the single scattered net flux at level L 
MS 

MS » Fnet : the multiple scattered net flux at level L 
Reldiff : relative fractional difference between the 
net fluxes (MS-SS) / SS 
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TABLE 6.7: 950 cm" 1 WINDOW BAND 
ZENITH ANGLES 60 # and 120* 


L 

SS„ „ 
Fnet 

MS_ 

Fnet 

Reldiff 

1 

0.2259E-04 

0.2310E-04 

0.0226 

2 

0.2026E-04 

0.2300E-04 

0.1352 

3 

0.1612E-04 

0.1976E-04 

0.2292 

4 

0.1335E-04 

0.1728E-04 

0.2944 

5 

0.1124E-04 

0.1543E-04 

0.3728 

6 

0.9550E-05 

0.1385E-04 

0.4503 

7 

0.8020E-05 

0.1244E-04 

0.5511 

8 

0.6530E-05 

0.5880E-05 

0.0902 

9 

0.5620E-05 

0.8050E-05 

0.4324 

10 

0.5280E-05 

0.8920E-05 

0.6894 

11 

0.4960E-05 

0.1011E-04 

1.0383 


L : level number 
SS 

SS ■ Fnet : the single scattered net flux at level L 
MS 

MS - Fnet : the multiple scattered net flux at level L 
Reldiff : relative fractional difference between the 
net fluxes (MS-SS) / SS 
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photons is so large as to overwhelm the downward directionality 
inherent in the aerosol particles which intercept the upwelling 
radiation. 

Another subtle physical process is at work which effects 
the number density of the aerosols as the zenith angle 
increases. The overall effect is for both MS and SS positive 
net fluxes (with MS positive fluxes still bigger) to decrease 
as the geometric path length is increased due to secant path 
elongation. The angular dependence is interpreted as effect of 
a secant stretching whereby the geometric path length is 
increased by 

Path length - s (l/cos9) (6.1) 

where (s) is the vertical distance between levels and (9) is 
the zenith angle. For example, the path length from the sur¬ 
face to 2.0 km looking at a 40® zenith angle becomes 2.6 km. 
Obviously, more scatterers will be included in the attenuation 
process to reduce the flux densities. But MS effects redirect 
photons back into the impinging pencil to increase over and 
above the SS process, which for LOWTRAN is a loss mechanism 
with no self-illumination. However, this process cannot con¬ 
tinue without abatement because of particle shadowing. 

The phenomenon of particle shadowing occurs when the 
number of intervening scatters along a path increases to such 
an extent that the physical separation between particles is 
greatly reduced from the perspective of traversing photon. Or 
put another way, the total volume of the scatters increases 
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within a column to the point where the scattered light from one 
particle is physically blocked by the presence of another. The 
opacity of such a path is greater and the flux densities would 
be reduced. Therefore, in the window region, increasing the 
number density of aerosol particles along a path will increase 
the diffuse component, when compared to first order scattering, 
but then it too is reduced as the particles get in the shadow 
of other adjacent particles. 






7.0 CONCLUSIONS 


7.1 Summary of Investigation 

The results of this investigation make it obvious that 
multiple scattering by aerosols of terrestrial infrared radia¬ 
tion is not without consequence. In juxtaposition over the same 
urban planetary boundary layer model (PBL) the multiply scat¬ 
tered total flux densities (MS), which involves scattering plus 
emission, are larger than singly scattered total flux densities 
(SS). This result is especially evident in the absorption band 
for water vapor (1200 cm ^ to 2000 cm *) where multiple 
scattering overwhelms single scattering by as much as 44 times 
(4,400%) in the lowest one third of the PBL. These vastly 
increased total flux densities are associated with water vapor 
and water coated aerosol emission emanating from the overlaying 
two thirds of the PBL where the visibilities are the lowest (< 
3.0 km) and the relative humidities are the highest (> 88.0%). 
In the water vapor window portion of the infrared spectrum (830 
cm ^ to 1250 cm ^), thermal emission from the surface is able 
to penetrate to the top of the urban PBL. Here, the relative 
fractional difference (MS - SS) / SS, is reduced but MS still 
is numerically superior ranging from 0.02 to 1.02 times (2% to 
102%) the SS flux densities. 
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In the window region at all levels the net flux densities 
remain positive (F ■ F+ - F+) for both single and multiple 
scattering processes. However, the distinction between the two 
diminishes with height as shown in Tables 6.5 to 6.7. This 
result occurs because particle shadowing is reduced allowing 
the multiple scattering process to dominate. This distinction 
in the lower levels is enhanced as the zenith angle is 
increased. This effect is caused by an increase in aerosol 
number density produced by secant path elongation. From the 
results it may be inferred that only at the lower levels (£ 
0.25 km) multiple scattering would enhance the positive flux 
densities and aid in the thermal emission from the surface. 

In the absorption bands of water vapor, negative total 
flux densities (F • “Ft - F+) predominate in the lowest one 
third of the PBL. Here, water vapor and water coated aerosol 
emission enhanced a downward directed atmospheric counter radi¬ 
ation as shown by the Tables 6.1 through 6.3. The appearance 
of negative flux densities below 0.50 km leads to the inference 
that absorption of photons is almost instantaneous and com¬ 
pletely masks any scattering process. It may, however, also be 
inferred that in those levels where water vapor absorption is 
less critical the scattering process dominates and the "cas¬ 
cade" of photons produced by numerous reabsorptions and reerais- 
sions from the layers above is released. The photons are 
directed downward by the aerosols that favor forward scatter¬ 
ing. When of an urban nature the aerosols have an asymmetry 
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factor of 0.85 and a backscatter fraction (8) of 0.15. It is 
concluded from these results that a synergistic relationship 
between absorption/emission and multiple scattering among the 
radiatively active gases and the aerosol is established, which 
produces an enhanced total downward directed flux that may aid 
in heating the lowest layers in an aerosol-laden PBL. 

7.2 Suggested Future Research 

The results of this investigation suggest that in the 
absorption bands of water vapor, in the presence of heavily 
aerosol-laden atmospheric structures, a greatly enhanced total 
downward directed net flux is produced. 

The anoriymous flux occurs at levels where the scattering pro¬ 
cess is able to dominate water vapor absorption. A recommenda¬ 
tion is made to study the possible influence this strengthened 
negative net flux would have on the temperature patterns within 
a PBL. This could be done by calculating the total flux diver¬ 
gence (AF . ) at levels located within the PBL and calculating 
net 

the radiative temperature changes by using the enthalpy rate 
equation presented by Liou (1980) as 

3T/3t - (-1/C p) (AF ./Az). (7.1) 

p net 

Infrared heating of this nature could play a role in a noctur¬ 
nal destablization of the lowest portion of the PBL. 
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